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SYMBOLS 
(aa  tlicy  appear) 

Instantaneous  velocities  in  the  x.  y,  and  z 
directions  respectively. 

Mean  velocities  in  the  x,  y,  and  z directions. 

Instantaneous  values  ;>f  the  velocity  fluctuations 
in  the  x,  y,  and  z directions 

Root-mean- square  values  of  the  velocity  fluctuation. 
Mean  pressure  at  any  point. 

Instantaneous  pressure  at  any  point. 

Instantaneous  value  of  the  pressure  fluctuation. 

Time. 


Air  density. 

Air  visroB  ity. 

Kinematic  viscosity  of  air. 

Kinetic  energy  of  the  turbulent  motion:  -jr  ft 

Instantaneous  value  of  the  speed  fluctuation 
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Local  half- angle  bet 
Reynolds  number: 
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Mean  Velocity  at  y = o. 

A pressure  term:  P + pv1 

Total  shearing  stress: 

‘X*uri8  the  value  at  the  wall. 


e 


uxr 


Distance  from  the  wall:  b - v 
Dimensionless  distance  from  *he  wall: 
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Skin  - friction  coefficient: 
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Fluid  weight  density. 

Height  of  fluid  column  in  manometer  - capillarity 
absent. 

Surface  tension,  force  per  unit  length. 
Dimensionless  meniscus  parameter. 

Additional  height  of  liquid  column  due  to  capillarity. 
Coefficients  in  King:s  Equation 
Hot-wire  Correction 


INTRODUCTION 


The  early  attempts  towards  a representation  of 
turbulent  shear  flow  were  by  the  30-cal!ed  phenomenological 
theories.  They  were  developed  during  the  late  twenties  and 
early  thirties,  and  /=*njr>yed  so me  success  at  nredictine  rnean- 
veiociiy  distributions,  but  it  became  apparent  that  their  basic 
assumptions  were  in  conflict  with  observations.  The  assump- 
tions and  discrepancies  were  clearly  reviewed  by  Eatchelcr 
(ref.  1). 


This  phase  of  development  emphasized  the  fact  that 
the  effects  of  turbulent  motion  (transport  of  stream  properties 
and  conversion  of  energy  to  heat)  required  the  investigation 
of  the  dynamics  of  the  fluctuating  motion.  The  simple  case 
of  isotropic  turbulence  had  been  studied  in  more  detail,  and 
some  success  has  been  realized  in  applying  the  results  of 
isotropic  theory  to  limited  portions  of  shearing  flow.  But  the 
production  of  turbulent  energy  from  the  mean  motion  and  the 
transport  of  mean  momentum  by  the  fluctuating  motion, 
important  to  shearing  flows,  is  absent  in  the  isotrop'c  case. 

In  fact  the  transport  of  all  mean  properties  may  be  different 
in  the  shearing  flow. 

It  became  apparent  some  years  ago  that  some  simple 
and  clearly  defined  turbulent  shearing  flows  should  be  set  up 
in  the  laboratory  in  order  to  accumulate  more  data  on  the  sta- 
tistical properties  of  the  turbulent  motion  preliminary  to 
formulating  a theory  of  the  mechanism  of  turbulent  motion. 

The  circular  jet  has  been  explored  by  Corrsin  and  Uberoi  (ref. 

2,  et.  al.  );  the  half-jet.  by  Liepmann  and  Laufer  (ref.  3);  the 
two-dimensional  wake  by  Townsend  (ref.  4,  et.  al.  );  the  channel 
(ref.  5)  and  pipe  (ref.  6)  by  Laufer;  the  boundary  layer  by 
Townsend  (ref.  7),  and  Klebanoff  and  Diehl  (ref.  8);  and  the 
boundary  layer  In  an  adverse  pressure  gradient  by  Schubauer 
and  Klebanoff  (ref.  9),  Ludwieg  and  Tillmann  (ref.  10),  and 
Clauser  (ref.  ll).  (The  last  two  were  restricted  to  measure- 
ments of  the  mean  motion  but  were  included  as  they  made  sig- 
nificant contributions  to  the  general  understanding  of  the  important 
parameters.  ) 


\ 
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The  boundary  layer  on  a solid  surface,  which 
is  important  in  engineering  and  convenient  to  explore 
experimentally,  is  characterized  by  a rising  Reynolds 
number  in  the  direction  of  the  flow.  This  implies,  as 
pointed  out  by  Townsend  < ref.  7 ) and  Clauser  ( ref.  11). 
that  a state  of  dynamic  equilibrium  cannot  exist  in  the 
turbulent  boundary  layer  (without  bleeding).  There 
would  be  a progressive  change  in  the  direction  of  the  flow 
between  the  large-scale  and  the  small-scale  structure 
of  the  motion,  e.  g.  the  ratio  of  an  "integral  scale"  to 
a "microscale".  This  complication  motivated  Laufer 
to  investigate  the  flow  in  a channel  with  a rectangular 
cross-eection  ( ref.  5 ) and  a pipe  ( ref.  6 ) where  the 
Reynolds  number  of  the  flow  far  from  the  entrance 
remains  constant.  The  latter  two  investigations  are 
closest  to  the  present  experiment,  so  will  serve  as  i’ne 
principal  basis  for  comparison. 


The  diverging  channel  is  a next  logical  investi- 
gation following  the  channel  and  pipe  since  the  character- 
istic overall  dimension  of  the  flow,  the  channel  width, 
is  varied  in  the  direction  of  the  flow,  but  the  dynamic 
equilibrium  is  retained  by  keeping  the  channel  essential" 
ly  two-dimensional  ( for  which  the  Reynolds  number  iB 
constant  ) and  the  walls  plane.  Reside  the  basic  sig- 
nificance of  such  an  investigation,  it  is  believed  that 
the  systematic  study  of  this  equilibrium  flow  affords 
the  soundest  approach  to  such  an  engineering  problem 
as  maintaining  unseparated  flow  in  a diffuser. 


Previous  measurements  had  been  reported  in  a 
diverging  channel  by  Donch  ( ref.  12  ) and  Nikuradse 
( ref.  13  ) but  they  were  limited  to  mean-velocities. 
Furthermore.  is  considerable  doubt  whether  the 

flow  in  reference  13  had  reached  equilibrium  ( discussed 
further  in  Appendix  I ).  Therefore  little  reference  will 
be  made  to  these  two  sources. 


The  two  most  important  forces  on  a turbulent  flow 
are  the  pressure  and  wall  shear  stress.  In  Laufer's  two 
experiments  these  were  in  balance.  Logically  one  might 
inve  f' ' gate  the  case  where  the  pressure  fovre  absent  — 
a very  small  positive-angle  diffuser  — but  it  was  ielt  that 
the  more  interesting  flow  is  that  where  the  pressure  forces 
are  large.  Therefore  a two-dimensional  diffuser  was 
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built  with  provision  for  a half-angle  of  divergence  from 
zero  to  6 degrees.  At  a 1 degree  half-angle,  which 
is  the  angle  reported  here,  the  pressure  force  was  7 
times  larger  than  the  wall-shear  force. 

Considerable  effort  was  spent  to  attain  an 
equilibrium  state  by  the  downstream  region  of  the  diffuser. 

A parallel-wall  channel  was  placed  ahead  of  the  diffuser 
to  give  the  flow  a head  start  towards  equilibrium,  and  the 
diffuser  test  section  is  117  half-widths  ( <^  , defined  later  ) 
long,  more  than  three  times  that  of  reference  13.  The 
steps  taken  to  establish  the  equilibrium  flow  are  discussed 
in  detail  in  Appendix  I.  The  flew  was  in  equilibrium, 
determined  by  the  mean  velocity  and  turbulence-level 
profiles,  in  a considerable  portion  of  the  diffuser. 

Townsend  ( ref.  7 ) first  demonstrated  the  import- 
ance of  the  turbulent  region  very  near  the  solid  surface. 

Later  Laufer  ronfirmed  the  observation  in  the  channel  (ref 
5 ) and  the  pine  ( ref.  6 ) where  he  demonstrated  that  the 
rate  of  production,  diffusion,  and  dissipation  of  turbulent 
energy  reach  a sharp  maximum  at  the  outer  edge  of  the 
laminar  sublayer  and  are  much  larger  than  in  the  other 
regions.  It  appears  that  this  region  deserves  intensive 
investigation.  This  contention  has  guided  the  work  in  the 
present  investigation.  The  diffuser  is  about  8 inches  wide 
at  the  downstream  end  and  very  small  probes  were  construct- 
ed to  enable  measurements  into  the  laminar  sublayer. 


ANALYTICAL  CONSIDERATIONS 


Equations  of  Motion  for  Flow  in  a Two-Dimensional 
Converging  or  Diverging  Channel  with  a Small  Angle. 

The  mean  - and  fluctuating  velocity  components  are 
represented  by  U(  and  Adi  respectively  in  X/  Cartesian 

coordinate  and  pressure  is  P.  Repeated  indices  indicate 
summation.  The  time  (or  ensemble)  average  is  indicated 
by  a bar,  so  for  stationary  flow 

Ufa,*)-  U;(X;)(-  UifriJ) 

pfc , i)  - p fa)  + pfc,  t) 

where  Alf  r pi  E 0 by  definition.  The  continuity 

equation  (1),  Reynolds  equations  (2),  and  turbulent  kinetic- 
energy  equation  (3)  are 


Vi. 
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t 

k. 
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Henceforth,  U,  V,  and  W will  replace  ^)‘  in  the 
respective  directions  x,  y>  and  z,  which  replace  TCi 
The  coordinate  x is  in  the  direction  of  flow  and  along  the 
line  of  symmetry;  y ia  the  transverse  coordinate  in  the 
ni2no  of  flow;  and  7 ifl  normal  thereto. 


Jr, 


For  flow  in  a two-dimensional  converging  or  di- 
verging (henceforth  "diverging"  will  be  used  to  imply  both) 
channel  at  a small  angle 

(a)  VV  *0 

(b)  C^L-^30  , where  ( ) is  any  term. 

d 2 

(c)  X I , where  S'-  is  the  local  half-angle 

between  the  walls. 

Equation  (1)  reduces  by  (b)  to 


iu 

dX 


*4) 
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The  following  assumptions  are  made  regarding  the 
turbulent  motion: 


and  u arc  u£  iiic  a aui 6 o*  u£x  of 


» 


f 

i 


’ 


(d)  U",  V,  W 
magnitude.  This  was  found  to  be  the  case,  except 
very  near  the  wall  where  if*  dominates,  and  at  the 
center  where  ti U"  goes  to  zero. 


— ; i it 

(e)  U'-  « Li 


inis  is  noi  universally  Vatiu 
for  example  in  the  outer  region 


of  a jet  entering  a fluid  at  rest,  near  a separation  point, 
or  possibly  in  a channel  at  sufficiently  low  Reynolds^ 
number.  However  Laufer  (ref.  5)  found  Ux  < O.l  \J 
everywhere  in  the  channel  at  Rl  > 3X/0+  - This 

was  confirmed  in  the  present  experiment. 

Since  Ck  « / it  is  reasonable  to  make  the  approxima- 

A. 

It 

d y 

where  £ refers  to  any  velocity  term.  This  will  be  called  the 
"boundary- layer  approximation".  Note  that  the  boundary-layer 
approximation  is  invalid  at  the  axis  of  symmetry.  In  this  region 
other  assumptions  will  be  called  upon.  Using  the  boundary- 
layer  approximation  in  equation  (4)  one  deduces  that  V « U . 

At  the  centerline  V«V  by  reason  of  symmetry,  so  it  applies 
to  the  whole  region  of  flow. 


The  x-momentum  and  y-momentum  equation  become 


tj  £11  + v = - 1 + y - 


6* 


d'i 


U A 


_£  A-  ( P + pir1)  -0 

P d<j  1 r ' 


(5a; 

(5b) 
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So  P +'  pV*  , called  p*  is  a function  of  x alone,  tty 

(H)  and  (e)  p{P)  is  small  in  equation  (5a)  so  it 

becomes  ^ * 

0 f V <tR  = - -L  d P*  , o d20  _ djuy 

dx  dy  f dx  d yl  d y 


The  turbulent  kinetic -energy  equation  (3)  by 
assumptions  (a),  (b),  (c),  and  (d)  is: 


I 

r 

$ 


1 
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The  fourth  term  is  negligible  compared  to  the  fifth  term  by  the 
boundary-layer  approximation,  except  at  the  channel  center 
where  the  fifth  term  goes  to  zero.  However  in  this  region  the 
fourth  term  is  negligible  compared  to  the  second  term  as  (e)  be- 
comes so  strong  that  U « U . Laufer  (ref.  6)  found  that  in  a pipe 
jj  was  important  very  near  the  wall  so  it  will  be  kept 

« J 

in  the  equation.  This  was  confirmed  in  the  present  case.  However 
even  there  it  does  not  dominate  the  equation,  so  the  term  y J e 

would  be  negligible  by  the  boundary-layer  approximation.  The 
latter  term  at  the  channel  center  is  small  compared  to  77 

Sx 

for  a high  Reynolds  number,  therefore  it  will  be  neglected.  * 

♦ The  term  ——  should  be  neglected  by  the  boundary- lay- 

o /v 


er  approximation.  This  was  substantiated  everywhere  with  the 
exception  of  a small  region  near  the  wall  where  (1  1(0'  > or 

4/*  , and  (2)  the  orientation  of  the  coordiantes  on  the  centerline 

make  s somewhat  larger  than  if  oriented  on  the  wall.  The 

latter  demonstrates  the  inherent  advantage  of  cylindrical-polar 
coordiantes  in  this  problem.  There  was  insufficient  time  when  this 
was  discovered  to  make  such  a revision,  and  furthermore,  the  term 
is  negligible  in  the  integral  form. 
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The  turbulent  kinetic-energy  equation  for  the  diverging  channel 
become  s 


AA  U 


u Z-Z  -h  V -4-  1 

dx  dy  dijl 

+yu/cL“;)f<L“i\ 

$y*  ' [dxjidxj 


(7) 


The  terms  are,  in  order:  the  production  of  turbulent  energy 
from  the  mean  motion,  convection,  diffusion  by  turbulent  motion, 
and  the  last  two  terms  combine  the  diffusion  by  viscous  forces 
and  dissipation  to  heat.  The  identification  of  the  last  two  terms 
is  given  in  reference  14.  Equations  (4),  (6),  and  (7)  comprise 
the  differential  equations  for  the  diverging  channel. 

Integral  Equations 

The  integral  equations  will  now  be  derived  after  von 
Karman's  method  for  mean  motion  (ref.  15  ) and  Cora  sin's  ex- 
tension to  turbulent  motion  (ref,  16  ).  Equation  (4)  upon  inte- 
gration across  the  channel  and  the  use  of  Leibnitz'  rule  becomes 
bfx) 

t-xpd!,-ubi±+K-V,  = o 

o 

where  b (x)  is  the  duct  half-width,  tib  aud  Vk  are  at  the  wall 
and  is  at  yaf?  . The  wall  no -slip  condition  and  symmetry 
of  the  flow  about  y *0  reduce  this  equation  to 

tfx) 

kjUdi  = 0 


(8) 


- 9 - 


By  the  same  reasoning  the  integral  of  the  x- 
momentum  equation  (6)  becomes 


i d - - _4  Tb 

VdxJUdy~  ptidx  + fUS 


(9) 


where  T 
energy  (7)  is 


Ihe  integral  of  the  turbulent 


the  only  terms  remaining  being:  the  production,  the 
longitudinal  gradient  of  turbulent  energy,  and  the  term 
entering  dissipation. 

At  this  point  an  interesting  comparison  can  be 
made  between  the  relative  magnitude  of  energy  going  to 
turbulence  and  the  energy  going  directly  from  the  mean 
motion  to  heat  for  the  flow  in  the  parallel  channel.  For 
this  purpose  the  so-called  mechanical  - energy  equation 
is  derived  by  multiplying  each  term  in  equation  (6)  by  fj 
and  then  integrating.  Following  the  above  procedure  one 
obtains : 


The  terms  are  respectively:  the  work  done  by  the  mean 
pressure  , the  longitudinal  decrease  in  the  flux  of  mean 
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energy,  the  energy  fed  to  turbulence,  and  the  direci  dissi- 
pation to  heat.  The  last  two-terms  are 

6(x ) i 

7 = fr  — i c)/7  C T.T.  . 

•r~  / * ■'***'/  r-=  CJ'V  - / • ~ — <3fT/ 

/ ~y  1 y“/»y  ti  < 

° 0 


where 


? = *.' 

>6 


rr-  -pM.v 


suppose  (1)  that  all  the  mean-flow  dissipation  and  turbulent 
energy  production  occur  within  a region  of  thickness  £ ad- 
jacent to  the  wall,  and  (2)  that  the  laminar  shear  stress  is 
approximately 


Zl_ 

T...  ' € 


= 0 


JL  ± ?}  ± l 
b • 


z-5  rL  * tv 

-0-y) 


in  the  <X.  = 0 channel. 


Evaluation  of  the  integrals  gives 


r - ILk 

~ 3 b 


2 

H 


X \ 

b 


This  means  that,  to  an  order  of  magnitude,  the  dissipation 
of  energy  directly  to  heat  is  as  large  as  the  energy  that  goes 
directly  to  turbulence. 


Longitudinal  Similarity 

The  integral  equations  will  now  be  further  reduced 
by  making  the  assumption  that  the  flow  in  succeeding  sections 
is  similar.  By  this  is  meant  that  distribution  of  all  the  time  - 
averaged  properties,  both  mean  and  fluctuating,  across  any 
section  is  identical  at  successive  sections  if  made  dimension- 
less by  the  appropriate  characteristic  quantities.  This  will 
be  called  "longitudinal  similarity",  and  it  is  in  contrast  to  the 
Karman  type  wherein  each  small  region  is  considered  to  be 
similar  to  every  other  small  region  even  if  the  regions  are  at 
the  same  longitudinal  station. 

The  argument  for  the  existence  of  simple  similarity 
is  the  following.  Consider  a region  in  the  plane  of  flow  bounded 
by  the  channel  walls  and  two  transverse  sections  a small  dis- 
tance apart  (an  isosceles  trapezoid).  The  equations  of  motion 
(Reynolds  equations)  are  characterized  by  a Reynolds  Number, 
which  can  be  formed  from  the  channel  half-width  , b.  and  spatial 
mean  longitudinal  velocity  across  the  section.  By  continuity, 
equation  (8),  this  Reynolds  number  has  the  same  value  for  all 
such  regions  in  a two-dimensional  channel  for  any  wall  configura 
tion.  Therefore  the  equation  is  identical  in  all  regions.  If  the 
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walls  are  plane  then  each  region  is  geometrically  similar  , 
that  is  the  physical  boundaries  of  each  region  are  also 
identical.  Therefore  one  expects  that  the  flow  would  be 
independent  of  the  particular  value  of?C  . There  are  two 
other  possibilities.  A longitudinal  periodicity  may  set  in, 
but  this  appears  unlikely;  or  the  flow  may  progressively 
change  until  it  reaches  some  state  of  breakdown  such  as 
separation.  The  latter  is  expected  to  happen  as  the  angle 
between  the  walls  is  increased. 


It  should  be  noted  that  the  flow  in  an  axially- symmetric 
diffuser  or  the  boundary  layer  on  a flat  plate  would  not  be  ex- 
pected to  exhibitlongitudinal  similarity  since  the  Reynolds  num  - 
bers  changes  in  the  direction  of  flow. 


It  is  also  worth  mentioning  that  if  the  laminar  sub- 
layer has  a characteristic  Reynolds  number  then  the  ratio  of 
the  thickness  of  the  aub-layer  to  the  half-width,  b,  in  the  two- 
dimensional  diffuser  would  be  constant  in  the  direction  of  flow, 
compatible  with  the  previous  remarks. 


Simple  similarity  implies. 


b 


where  we  recall  The  following  integrals  will  be 

useful: 


I,* ft  It - f f'M  h = JffdV 


It  -jt fhd’l  l S -jhai^i 


7*  „ r . / / 

u-j  jnciy 

o 
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The  equations  of  continuity  (8),  momentum  (9),  and  turbulent 
energy  (10)  become  respectively 


lt  = A CONSTANT 


(12) 


/)//* 
/ vo 


=Z  JZ0C  + 


PUJ 


T r~.  , 

x5  I b dec 
RbL2  c/p< 


where  <K  is  the  diffuser  half-angle.  Let  us  proceed  without 
making  the  restriction  of  geometric  similarity,  i.  e. 

SjJ?  - A CONSTANT  , to  see  whether  similarity  could  exist  for 

a more  general  wall  shape.  Generally 

If  there  are  no  large  gradients  of  OC  then  for  a large  Rey- 
nolds number  the  turbulent  energy  equation  becomes 


(14) 


One  can  now  determine  the  conditions  for  longitudinal  similarity 
by  two  methods. 


(1)  Consider  equation  (13).  l or  smooth  surfaces 
and  similai  profiles  — . would  be  a constant.  Clauser 

pu* 

y fa  length)  of  V* 
Tu,  dX 


(ref. 


li 


presents  arguments  to  show 


would  have  a constant  value  for  similar  profiles  in  the 
turbulent  boundary  layer.  In  the  present  case  all  lengths 
would  remain  proportional,  therefore  b is  selected.  From 
this  equation  (13)  shows  would  be  a constant,  that  is 

the  walls  must  be  plane. 


i 

f 
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EQUIPMENT 


Wind  Tunnel 

A blowing-type  wind  tunnel  was  constructed  for  the 
present  investigation  ( figure  1 ).  The  room  air  passes 
through  a rloth  filter,  3 axial  flow  fans,  and  a settling 
chamber.  This  and  the  remainder  of  the  tunnel  ia  67in.  high. 
The  stream  then  passes  through  a lin.  wide  parallel-wall 
channel  and  enters  the  test  section.  The  fan  speed  was 
controlled  to  keep  a consta.nt  Reynolds  number  in  the  test 
section  with  day-to-day  changes  in  density  and  viscosity  of 
the  air. 


The  test  section  is  loft.  long.  The  side  wails  swing 
laterally  on  a sheet-metal  hinge  attached  to  the  parallel  section. 
Considerable  care  was  taken  to  make  the  walls  smooth  and 
plane.  Each  w<all  is  a symmetrical-sandwich  type  construction 
with  1 /4in.  tempered  masonite  faces  and  3/4in.  marine  ply-- 
wood  ribbing  on  lit.  spacing.  Surveys  were  made  of  the 
channel  width  at  = O over  a wide  range  of  atmospheric  con- 
ditions { temperature"  and  humidity  ) and  the  variations  in  the 
lengthwise  direction  were  less  than  1:500  ( width  variation/ 
unit  length  ) which  is  equivalent  to  0.  1 degree.  The  pressure 
drop  was  measured  at  =0  and  b = J.  in.  The  computed 
value  of  cf  was  20%  higher  than  I.aufer  (ref.  5 ) reported,  but 
only  6%  higher  than  Skinner  ( ref . 17  ) measured  later  in  the 
same  tunnel  as  ref.  5.  Since  a 2%  error  in  b would  account 
for  the  latter  difference  it  was  not  pursued  further.  The  walls 
were  then  waxed  and  polished,  and  are  dusted  before  each  run. 

Considerable  effort  was  spent  attaining  acceptable  two- 
dimensionality  of  the  flow  in  the  test  section.  The  channel  is 
pressurized  by  a throttle  at  the  outlet,  and  air  is  bled  ( an  esti- 
mated one  percent  of  the  total  flow  ) from  all  four  corners  by 
0.G10  in.  gaps  running  the  full  length  of  the  test  section.  In 
addition  a laminar  boundary-layer  tripper  was  installed  in  the 
parallel  section  to  remove  more  momentum  from  the  mid-plane 
by  hastening  the  onset  of  turbulence.  The  whole  program  is 
discussed  in  detail  in  Appendix  I.  But  the  final  result  is  that  even 
with  the  above  steps  the  growth  of  the  floor  and  ceiling  boundary 
layer  causes  some  momentum  to  be  convected  intothe  mid-plane 
by  vertical  flow  so  that  the  two-dimensional  momentum  equation 
is  in  error  by  about  10%  of  the  largest  terms  at  the  furthest 
downstream  station.  Total  pressure  probes  were  installed  per- 
manently at  three  elevations  of  Lite  furthest  downstream  station 
and  the  flow  is  regularly  checked  during  each  run.  A 5ft.  wall 
was  added  to  the  downstream  end  of  the  test  section  to  extend  the 
diffuser  and  reduce  upstream  effects  of  the  throttle. 
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Traversing  Apparatus 

The  probes  arc  posited  by  a traversing  mechanism 
with  a micrometer  head  graduated  in  0.  CGI  in.  , similar  to  the 
one  described  in  reference  5.  The  distance  between  the  probe 
and  its  image  reflected  in  the  w?11  wa  a measured  by  a tele- 
scope during  each  run  so  that  the  distance  from  the  wall  was 
estimated  tc  be  accurate  to  0.  001  in.  for  the  impact  probe  and 
0.002  in.  for  the  single-wire  hot  wire.  The  hot-wire  holders 
could  be  rotated  in  the  horizontal  and  vertical  planes  to  align 
with  the  stream  direction.  The  line  of  traverse  was  always 
normal  to  the  centerline  of  the  channel. 

Mean  Velocity  and  Pressure  Instruments 

A flattened  hypodermic  needle  was  used  for  most  of 
the  mean -velocity  measurements.  The  opening  was  0.  007  in. 
x'0.042  in.  and  the  wall  thickness  was  0.  004  in.  Near  the 
surface  several  measurements  were  made  with  a 0.  0005  in. 
dia.  platinum  hot-wire.  The  wire  was  5/16  in.  long  anu  held 
taut  by  two  small-diameter  needles.  The  mean  resistance  of 
the  wire  was  kept  constant  throughout  the  flow  traverse. 

Static  pressure  was  measured  by  static  taps,  each 
made  of  a l/’8in.  dia.  rod  threaded  into  the  waii  and  held  securely 
by  a jam  nut  glued  into  place.  The  orfice  was  0.  020  in.  in  dia. 
and  finished  smooth  on  the  air  flow  side.  There  are  51  static 
taps  and  7 traversing  points  in  all. 

Due  to  the  low  velocities  in  the  downstream  portion  of 
the  test  section  it  was  necessary  to  develop  a very  accurate 
pressure  measuring  instrument.  This  is  discuseed  in  detail  in 
Appendix  II.  Water  and  a wetting  agent  were  used  in  an  inclined 
manometer  with  an  18:1  slope.  By  locating  the  meniscus  with  a 
traveling  microscope  an  accuracy  in  pressure  of  10”^  in.  of 
water  was  attained,  which  is  quite  satisfactory.  An  attempt  was 
made  to  reduce  the  meniscus  motion  due  to  pressure  fluctuations 
by  placing  constrictions  in  the  connecting  lines  so  the  whole 
system  had  a time  constant  of  1 minute.  The  fluctuations  of  the 
meniecua  were  sufficiently  damped  to  identify  mean  values  with- 
in the  desired  accuracy,  but  the  effect  of  the  non-zero  influx 
and  efflux  at  the  probe  is  still  in  question. 

Hot-Wire  Equipment 

The  compensating  amplifier  and  associated  control  unit, 
calibration  unit,  and  power  supplies  arc  essentially  duplicates 
of  those  described  by  Kovasznay  in  reference  18.  It  is  a constant- 
current  type.  Compensation  for  the  hot-wire  thermal  lag  is 
adjusted  manually  using  the  square-wave  technique.  In  the  present 
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case  a single  amplifier  channel  was  employed.  Of  importance 
to  this  investigation  is  the  fact  that  the  amplifier  frequency 
range  extends  down  to  2 cycles  per  second.  The  amplifier  out- 
put filter  was  nn  position  "A"  which  has  a response  of  1 3 

percent  from  2 to  2500  cycles  per  second.  At  station  5 the 
signal  was  submerged  into  the  noise  between  2000  and  3000  cps. 

The  calibration  of  the  amplifier  is  shown  in  figure  3. 

All  output  readings  were  taken  by  a vacuum  thermocouple 
and  millivoltmeter.  A block  diagram  is  shown  in  figure  4.  The 
Lr'  , US1 , and  Utr  signals  were  fed  directly  from  the  amplifier 
to  ihe  vacuum  thermocouple. 

An  R-C  circuit  was  found  most  suitable  for  differentiation. 

The  complete  differentiating  circuit  was  linear  from  20  to  1000 
cycles  per  second  and  had  only  a 15  percent  departure  at  2000 
cycles  where  the  hot-wire  signal  generally  was  at  the  noise  level. 

Spe  cirum  measurements  were  made  using  a Hewlett-Packard 
Model  300A  harmonic  analyzer  with  a nominal  frequency  range  of 
30  to  16,  000  cycles  per  second.  A cathode  follower  was  added  to 
the  output  stage,  to  accommodate  the  relatively  low  resistance  of 
the  vacuum  thermocouple.  All  measurements  were  made  with 
the  30  and  145  nominal  half-band  widths. 

Platinum  wires  made  by  the  Wollaston  process  were  UBed 
exclusively  for  the  probes..  They  were  etched  then  soft  soldered  to  the 
prongs.  For  u'  in  the  channel  center  the  wires  were  nominally 
0.  00015  in.  in  dia.  , 0.  04  in.  long,  and  attached  to  fine  sewing 
needles;  near  the  wall  the  wires  were  0.  01  in.  long  and  attached  to 
fine  jeweler's  broaches  to  reduce  aerodynamic  interference.  Also 
to  measure  near  the  wall  the  x -me ter  probe  holder  used  for  all 
the  U*  ,U//,  and  ly-  measurements  was  constructed  of  fine  jewel- 
ers broaches.  They  were  spaced  0.  025  in.  apart  and  0.  015  in. 
between  pairs.  Each  wire  was  0.00015  in.  in  dia.  and  inclined 
approximately  40  degrees  to  the  free  stream  to  increase  the  later- 
al sensitivity  compared  to  the  longitudinal.  The  time  constant 
for  all  the  above  probes  was  0.  5 to  0.  8 milliseconds. 

For  measurement  of  the  a 0.  00005  in.  dia.  wire 

dt 

was  used.  Its  time  constant  was  0.1  milliseconds  and  the  noise 
was  reduced  to  2 percent  of  the  output  signal,  compared  to  10 
percent  when  differentiating  IT  and  ur  as  measured  with  the 
larger  wires.  The  same  wire  was  used  to  measure  the  spectrum 
nearest  the  wall.  All  wires  were  "overheated"  50  percent  above 
the  cold  resistance. 
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Preliminary  Measurements  in  1 in.  Channel 

The  preliminary  measurements  in  the  1 inch  channel  were 
made  at  9L  - 0,  The  longitudinal  static-pressure  distribution, 

mean-velocity  distribution,  and  centerline  trubuler.ee  level  were 
measured  at  = 12,200  and  compared  to  Laufer^s  values 

(ref.  5)  in  the  straight  channel. 


The  static-pressure  distribution  is  shown  in  figure  5.  The 
gradient  is  quite  smooth.  The  computed  value  of  was 

somewhat  higher  than  L*»ufer  obtained  but  the  difference  is  probably 
within  the  limits  of  the  experimental  accuracy  (see  p.  19). 


The  mean-velocity  distribution  at  - 2 83  is  compared 

in  figure  6 to  L,aufer  s distribution.  The  values  nearer  the  wall  tend 
to  be  high,  but  this  may  be  attributable  to  the  difference  between 
the  hot  wire  (ref.  5)  and  the  total-head  probe  (present  experiment), 
although  the  correction  of  page  21  applied  to  the  hot-wire  measure- 
ments of  reference  5 increased  jj  at  = o.l  by  only  C.6  per- 

cent. The  total-head  distribution  was  repeated  several  times,  but 
the  profiles  duplicated.  The  cause  of  the  difference  is  an  unsettled 
question . 

aj.'  "i  ’ 

The  turbulence  level  was  -g-  - .036  at  -jr  = 1,  com- 
pared to  = .035  for  reference  5.  The  difference  is  within 

the  experimental  accuracy. 
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PROCEDURE  AND  RESULTS 

Mean  Velocity  and  Static  Pressure 

All  measurements  were  mado  at  U.  = 1 deg.  and 
= 31,  400,  the  latter  value  to  compare  to  one  of  Laufer'a 
experiments  at  06  = O deg.  ia  125  feet  per  second 

at  the  inlet  c i the  test  section  and  decreases  to  16  feet  per 
second  at  station  5. 


The  mean  static-pressure  distribution  along  the  wall 
of  the  diffuser  at  the  midplane  is  shown  in  figure  7.  The  slope 
of  the  faired  line  corresponds  to  = 0.  010. 

‘tra* 


The  pressure  was  measured  by  the  difference  between  stations 
along  the  test  section. 


The  mean-velocity  distribution  is  shown  in  figures  a and 
9.  At  stations  2,  3,  and  4 traverses  were  made  in  the  center 
of  the  channel  for  similarity  comparisons.  The  most  complete 
traverse  was  made  at  station  5 where  all  of  the  data  near  the 
wall,  figure  9,  were  obtained.  The  nearly  - straight  line  from 
the  wall  in  figure  9^-^pre sents  the  (laminar)  velocity  distribution 
computed  from  <lsL  — <)P  (the  inertia  term  is 

r<hym  Sx  negligible  in  this  region) 

where  at  the  wall  corresponds  to  Tu  = 0.0013, 

arrived  at  later  in  this  report. 


The  profile  was  symmetrical  on  opposite  sides  of  the 
channel  at  station  5 within  1/2  percent  of  •M.  The  accuracy 

J7 

in  measuring  y decreases  downstream.  At  station  5 the 
probable  error  f 0.  675  x std,  deviation  ) of  (7  is  estimated 

- 0.  1 diminishing  towards  the  center. 


to  be  0.  005  at  ^ 

h 


Most  of  the  mean  velocities  were  computed  from  the  total- 
head  probe  measurement  and  the  wall  static  pressure  to  which 
the  following  corrections  were  applied; 


1.  ) Non-linear  response  of  probe  to  magnitude  and 
direction  of  velocity.  The  total-head  probe  was  tested  for 
directional  sensitivity,  and  with  the  probe  axis  oriented  in  the 
x-direction  it  was  found  essentially  insensitive  to'/  and  W over 
a range  of  - 30  deg.  in  each  plane,  (i.  e.  the  probe  has  a cosine- 
squared  response  to  rotation  in  a uniform  stream.  ) 
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To  estimate  the  effect  of  velocity  fluctuations  we 
will  limit  our  considerations  to  the  external  flow  at  the  probe, 
i.e.  , we  will  ignore  the  dynamics  of  the  f'uid  within  the  probe, 
tubes,  and  manometer.  We  will  also  assume  that  the  pressure 
response  at  the  probe,  (pt('t)  : is  in  equilibrium  with  the 
local  instantaneous  velocity  and  pressure.*  Then  neglecting 
size  affects  and  including  the  above  directional  response 


R(t)  = ieu'(t)  * Pit ) 


= 4-pfu 
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Define  an  apparent-mean  velocity,  \j 


• by 
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Fbr  a low  turbulence  level 
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was  measured 


which  was  the  form  employed,  where  tti 
by  the  hot-wire  anemometer.  The  correction  was  in  the  order 
of  5 percent  or  less. 

Z.  ) Finite  probe-  width.  Essentially  the  correction 
of  Young  and  Maas  ( ref.  19  ) was  employed  to  account  for  the 
-finite  probe  width  in  a transverse  velocity  gradient,  except 
that  the  minimum  probe  width  was  used  in  place  of  the  tube 
diameter.  The  correction  takes  the  form  of  a shift  in  the 
apparent  centerline  of  the  probe.  The  maximum  shift  was 

=a  C,  00065", 

is  which  is  a large  correction  in 

the  laminar  sublayer.  The  presence  of  the  wall  exposes  this 
correction  to  some  question. 

3.  ) Reynolds  number  of  prooe.  Due  to  the  very  low 
velocities  near  the  wall  the  viscous  effects  around  the  probe 
were  examined.  A correction  was  made  following  the  method 
of  Hurd,  Chesky,  and  Shapiro  { ref.  20  ),  where,  again,  the 
minimum  width  of  the  probe  was  employed  in  place  of  a tube 
diameter.  The  maximum  correction  was  2 percent. 


♦ The  quasi-equilibrium  analysis  requires  justification. 
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Equation  (5b)  indicates  a transverse  difference  in  the 
pressure  due  to  . This  correction  to  the  mean  velocity 

would  be  nearly  everywhere  in  the  order  of  0.  6 percent  going 
to  zero  towards  the  wall.  It  was  considered  negligible  compared  to 

M-  . 

I 1 

In  view  of  the  large  corrections  to  the  total-head 
measurement  near  the  wall  the  mean  velocity  was  also 
measured  with  a hot-wire.  The  nonlinearity  of  the  latter  is 
such  that  the  values  tend  to  be  low,  so  the  following  correction 
was  employed: 


The  method  will  be  merely  outlined  since  it  is  considered 
tentative  pending  an  analysis  for  large  fluctuations  as  existed 
in  the  laminar  sublayer.  The  method  includes  the  other  leg 
of  the  bridge  and  essentially  fits  a parabola  to  King’s  law 
(current  vs.  velocity)  about  the  point  of  the  mean  velocity  for 
the  response  to  fluctuations.  This  assumption  is  limited  to 
small  fluctuations.  For  a wire  oriented  in  the  z direction  the 
"apparent"  mean  velocity  (corresponding  to  the  mean  current 
and  mean  resistance  in  King's  equation)  is 


/ + 


Ul 


I O' 


where  A and  B are  the  constants  in  King's  equation 


I*R 

f?  - K 


The  wire  hot  and  cold  resistances  are  R and  JRA  respectively. 
The  correction  was  in  the  order  of  five  percent  or  less. 

It  is  to  be  noted  that  the  velocities  measured  by  the  hot 
wire  also  fall  below  the  curve  obtained  from  the  wall  shear 
stress.  The  wire  was  calibrated  at  velocities  in  excess  of  10 
feet  per  second  and  used  down  to  1 foot  per  second,  so  an 
error  would  have  been  incurred  in  extrapolation.  Also  a 
possibility  of  error  exists  due  to  the  determination  of  y‘ - o; 
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0.003  in.  would  put  the  points  cn  the  curve.  The  difference 
between  the  wire  positions  could  be  measured  very  accurately, 
however,  so  the  velocity  gradient  was  relied  upon.  Moreover, 
two  extreme  extrapolations  ot  the  hul-wiie  calibration  changed 
the  computed  gradient  by  only  two  percent. 

In  a recent  private  communication  Mr.  Richard  Cox 
indicated  that  some  measurements  he  made  in  Gottingen  show: 

(a)  that  the  correction  for  velocity  fluctuations  is 
larger  than  computed  above, 

{b}  the  conduction  to  the  wall  is  measurable  to  at  least 
50  wire  diameters  from  the  wall  (y’/b  = 0.  006). 

The  exact  magnitude  of  these  corrections  is  not  known  more 
definately  by  this  writer  so  they  could  not  be  applied.  However, 
qualitatively  they  do  account  for  the  difference  of  the  measured 
points  from  the  curve.  It  is  planned  to  investigate  this  further. 

Finally,  King's  analysis  also  indicates  that  the  hot-wire 
equation  changes  form  at  a velocity  equivalent  to  U - 0.1, 

a 

where  the  measured  values  change  slope. 

The  above  comments  focus  attention  on  the  importance  of  an 
intensive  investigation  of  the  instrument  calibrations  at  low 
velocities,  near  a surface,  and  with  a large  gradient  in  the  mean 
velocity  in  view  of  the  importance  of  the  laminar  sublayer  to 
the  whole  flow. 

The  establishment  of  a fully-developed  flow  was  a prime 
objective  of  this  investigation.  A comparison  of  the  mean- 
velocity  profiles  in  figure  6 indicates  that  some  peaking  of  the 
profile  did  take  place  downstream  irom  station  2.  But  because 
of  the  considerable  length  of  the  test  section  it  was  possible  to 
determine  that  the  change  in  the  profile  between  stations  2 and 
5 was  insignificant.  This  is  presented  in  detail  in  Appendix  I. 

It  is  evident  from  figure  8 that  the  difference  in  the  mean 
velocity  destribution  between  o(.  = 0 deg.  and  oC  = i deg. 
is  substantial.  This  is  not  surprising  inasmuch  as  the  force 
due  to  the  longitudinal  pressure  gradient  has  increased  to  a value 
7 times  larger  than  the  wall-shear  stress  force  (which  has  changed 
only  slightly). 
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Turbulence  Levels  and  Shearing  Stresses 

The  single-wire  probes  for  measuring  U*  were 
calibrated  by  the  constant  resistance  method  of  Kovasznay 
(ref.  21)  and  used  at  essentially  constant  current.  Laic*  U* 
was  checked  by  the  x-meter  which  was  calibrated  with  the 
current  constant,  as  used,  and  the  values  were  found  io  be 
25  percent  higher  over  most  of  the  channel.  The  trouble 
was  finally  traced  to  the  dependence  of  the  so-called  constant, 

A,  in  King’s  equation  ( ref.  22  ) on  the  wire  temperature. 

It  is  attributable  to  the  variation  in  the  thermal  conductivity 
of  the  air  surrounding  the  wire  with  th*»  wire  temoerature. 

This  effect  had  been  identified  and  experimentally  evaluated 
by  King,  but  evidently  has  been  overlooked  in  a number  of 
intervening  turbulence  investigations.  Since  everyone  who 
measures  s»i£ficiently  close  to  a surface  is  faced  with  this 
effect,,  as  it  incresses  at  low  velocities,  it  will  be  presented 
here.  The  full  derivation  is  given  in  Appendix  III,  and  will 
be  merely  outlined  below. 

King’s  equation  can  be  written 

IXR  - {a(* ) * BVu](-»  - i9e) 

where  I*  R is  the  heat  input  to  the  wire  in  equilibrium; 
is  the  so-called  "constant",  which,  for  a given  wire,  exhibits 
the  dependence  upon  the  wire  temperature  is  the 

air  temperature,  assumed  constant;  B,  the  velocity  coefficient, 
assumed  constant,  is  only  slightly  temperature  sensitive. 

King's  empirical  expression  for  A <*>  is 

A - A*  [ / -+  - T°e)7 

where  C = 0.  00114  deg.  cent7* 

Inc’uding  the  second-order  coefficient  in  the  thermal 
resistivity 

* ” R,  [/  * <(■&  -■»%)■*  &(■&  - 

where  ia  the  reference  temperature  for  R,  and  for  constant - 

current,  I,  operation 
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where,  very  nearly, 
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and  the  ’’overheats"  are  defined 
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in  which  Re  ami  Rf  correspond  to  and  "0^*  • For  the 

case  where  A is  a constant  and  ^ = © , J*  - O , and  the 
equation  of  ref  21  is  obtained.  The  value  of  was  0.  16  to 
0.  24,  increasing  towards  the  wall.  After  this  correction  was 
applied  the  single-wire  values  increased  to  within  5 to  10  percent 
to  the  x-meter  values,  the  agreement  being  better  near  the  wall. 
In  a^private  communication  Lauler  indicated  that  he  had  used 
CL  hr'  = 0,  2 in  the  channel  experiment  and  the  velocities  were 
somewhat  larger  at  a corresponding  Reynolds  number,  so  the 
corrections  for  his  data  would  be  in  the  order  of  5 percent.  This 
correction  was  not  applied  to  his  data. 

The  values  of  IT  are  given  >n  figure  10  and  11 . The 
profiles  at  stations  2 and  5 are  compared  in  figure  10  in  order  to 
verify  that  the  turbulent  motion  wa.s  also  fully-developed,  as  was 
the  mean  velocity.  It  is  interesting  that  the  maximum  value 
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near  fhe  wall  is  nearly  equal  to  Laufer's  maximum  value,  in 
particular  if  the  above  corrections  were  applies,  although  some- 
what displaced,  indicating  perhaps  that  the  maximum  value  of 

u 

is  more  sensitive  to  the  Reynolds  number  than  to  the  pressure 
gradient. 

The  method  of  measuring  IT  , t» T*  , and  WV*  is 
similar  to  the  procedure  described  in  reference  23.  The  probe 
was  placed  in  a uniform  stream  and  rotated  about  the  neutral 
position  to  determine  the  «/*  - sensitivity  ( or  u/*  - sensitivity  ), 
and  then  the  velocity  was  varied  to  determine  the  U -sensitiv- 
ity, all  at  essentially  constant  current.  This  method  is  con- 
siderably less  laborious  than  calibrating  the  wires  separately, 
and  introduces  less  possibility  of  error.  Also  the  probe  is 
cydibrated  in  the  circuit  exactly  as  used.  With  two  identical 
wires  the  ■'neutral  position",  at  which  the  surn  and  difference 
of  the  voltages  across  the  wires  would  correspond  simply  to  u 
and  L r*  fluctuations  respectively,  could  be  easily  determined. 
Actually,  since  the  wires  inevitably  differ  in  diameter,  length, 
and  straightness  it  is  necessary  to  make  compensations.  This 
is  accomplished  as  follows.  Consider  the  voltage  sensitivity  of 
wires  no.  1 and  no.  2 to  xiuctuatious  11  and  U"  to  be  of  the 
form 


e» 

- a,  u 

. 1 . ^ 

+ O/tT 

e* 

=■  ^2  M 

W bz  KT 

Then  if  a^  = 

and  b^  = b? 
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Actually  three  conditions  were  satisfied,  the  third  being  that 
the  wires  had  the  same  overheat.  We  have  3 degrees  of 
freedom  available  to  satisfy  these  conditions,  provided  the 
differences  are  small,  namely: 

(1)  the  "tap-off"  voltage,  i.e.  a potentiometer  to  reduce  trie 
measured  voltage  from  one  wire,  (2)  the  current  into  each 
wire  by  adjustment  of  the  resistance  in  series  with  the  wire, 
and  (3)  the  angle  of  the  probe  to  the  air  stream.  It  turned  out 
that  for  the  probe  used  only  adjustments  (1)  and  (3)  were  necessary 
to  obtain  overheating  values  within  6 percent. 

The  procedure  was  to  approximate  the  neutral  angle 
and  measure  aj  - a 2 keeping  the  direction  constant,  and 
bj  - hz  keeping  the  velocity  constant.  The  latter  was  the 
most  difficult  to  measure  with  certainty  and  consequently  set 
the  limit  of  the  final  accuracy;  the  difficulty  was  in  maintaining 
the  velocity  sufficiently  constant  that  ej  + ez  changed  due  only  to 
IT  . By  knowing  the  effects  of  changing  the  probe  angle,  etc.  , 
which  were  put  into  an  equation,  it  was  possible  to  change  the 
angle  and  tap-off  voltage  simultaneously  so  that  the  conditions 
were  satisfied. 

— “*"2 

In  order  to  measure  UW  and  KJ  with  a 5 percent 
accuracy  the  a-coefficients  must  be  matched  to  approximately 
2 percent  and  the  b-coefficients  to  4 percent.  This  corresponds 
to  setting  tile  probe  angle  to  about  1 degree.  It  is  believed  the 
neutral  angle  could  be  determined  to  i.  1-2  deg. 

King's  law  was  used  to  extend  the  directional  sensitiv- 
ities to  the  lower  mean  velocities,  but  no  $ -correction  was 
necessary.  The  Ut7  signal  had  large  low  frequency  fluctuations, 
since  it  is  in  the  non-isotropic  portion  of  the  spectrum,  which 
added  to  the  difficulty  in  determining  the  mean  values  for  Uu- 


but.  CT*  readings  were  quite  3table. 

, • 

j j hf 

The  distribution  of  — - and  - are  given  in 

■ — u u 

figure  10,  and  --SC.  in  figure  12. 

uu'  

The  absolute  values  of  UlT  are  shown  in  figure  13.  The 
curve  for  "C  was  computed  from  the  pressure  gradient  and 
mean-velocity  distribution  using  a form  of  equation  (b)  wherin 
longitudinal  similarity  was  assumed.  {A  correction  was  applied 
due  to  the  vertical  flow  of  momentum  by  adjusting  the  value  of 
pressure  gradient  by  10  percent  so  Tur  = 0.  0013,  the 

P 

value  arrived  at  below.  ) The  dashed  lire  was  obtained  by  sub- 
tracting the  viscous  shear  stress  and  pressure  gradient  terms. 
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Of  importance  to  later  correlationa  is  the  wall  shear  stress. 
This  was  computed  from  three  independent  measurements: 

1.  x -meter.  Using  the  best-fitting  curve  of  UtT  at 

9' 

i =0.1.  figure  13,  the  value  at  the  wall  was  computed  by  the 


= 0.  0014. 


momentum  equation  to  be  = 0.  0014. 

cW- 

2.  Hot-wire.  The  mean-velocity  gradient,  figure  9, 
measured  by  the  the  hot  wire  in  the  outer  region  of  the  laminar 
sub-layer,  corrected  for  the  difference  ( 4 percent  ) from 

due  to  pressure  gradient,  gave  = 0.00133 

<o  4* 

Possibly  this  is  less  accurate  tharTl  in  view  of  the  effect? 
mentioned. 

3.  Total-head  probe.  This  is  the  most  questionable 

since  the  probe  is  0.015  in.  wide  compared  to  a laminar  sub- 
layer of  about  0.  04  in.  The  best  laminar  profile  through  the 
three  inner  points  gives  *£*«*?  = .0011. 

c* 

_ _ The  final  value  taken,  the  average  of  the  three,  was 

= 0.0013.  This  is  estimated  to  be  accurate  to  10  percent. 


The  absolute  turbulence  values  are  compared  in  figures 
14  and  15.  Important  is  the  fact  that  the  turbulent-energy  level 
is  roughly  4 times  larger  than  the  cC  = 0°  case  over  a large 
region  in  the  center  of  the  channel,  decreasing  to  become  equal 
near  the  wall.  This  is  a substantially  greater  effect  of  the 
pressure  than  was  exhibited  on  the  mean  velocity.  Also  it  is 
evident  that  the  larjge  scale  motion  is  not  isotropic  as  l/*  is 


roughly  equal  to  t/* 


in  the  central  region. 


It  can  be  shown  from  the  turbulent  energy  equations  £ar 
the  individual  components  that_the  production  occurs  in  the  * 
equation  and  energy  is  fed  to  If2  and  U/’a  through  the  static 
pressure  fluctuations.  Figure  14  shows  that  the  energy  is 
divided  equally  between  If*  and  Uf * in  the  center;  towards 
the  wall  IT*  is  smaller  due  to  the  presence  of  the  surface. 
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Spectrum 

The  une  dimensional  energy  spectra  of  U , 1/  and 
Uf  are  shown  in  figures  16,  17,  and  18.  The  signals  were 
calibrated  by  comparing  the  output  at  several  frequencies  in 
each  profile  to  that  for  a known  sine-wave  input.  An  "effective" 


(rectangular } banu  width  uf  13  i..  y,  a.  was  used. 


mL  • 1..  .« 

X mu  VC^AV4V> 


was  obtained  by  measuring  the  band-pass  shape,  and  was  con- 
flrmed_by  comparing  the  u*  by  integrating  the  whole  spectrum 
to  the  u*  value  measured  directly.  This  was  carried  out  in 
another  wind  tunnel  at  a velocity  high  enough  that  the  energy  at 
lower  wave  numbers  could  be  measured.  Comparing  the  band- 
passed  signal  for  a white-noise  input  would  be  more  dependable, 
but  a suitable  noise  source  was  not  available  at  that  time. 
Therefore  the  absolute  level  of  the  energy  spectrum  is  somewhat 
in  question,  but  the  shapes  of  the  spectra  should  be  reliable. 

7fhe  t(  - spectra  were  corrected  for  the  wire-length  effect  by 
the  method  of  reference  24. 


the  fraction  of 


The  spectrum  is  defined,  for  example, 
t?  at  wave  number 


by  i-Jfc.)  being 
and  normalized  by 


The  microscales  were  computed  from  the  spectra  by 
the  well-known  relations 
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They  compare  very  well  with  the  values*  obtained  by  differenti- 
ation {described  below),  figure  21. 


Production  Terms 

In  the  Cartesian  coordinate  system,  as  mentioned  in 
the  Preliminary  Analysis,  the  following  are  two,  of  the  three 
production  terms  { the  remaining  term ,_«77J,  J2M  ie  ne8“ 
ligible  ) fa 

_ jrr.  / “5 

- Uu-  - (u-  ~ 


The  second  term  is  small  everywhere,  except  near  the  wall, 
if  the  coordinates  are  parallel  to  the  centerline  as  defined  in 


this  investigation.  If  the  coordinates  are  aligned  parallel  to 
the  wall  it  is  small  even  there,  (Of  course,  the  tirst  tcxul  i*t“ 
creases  simultaneously  as  the  sum  remains  constant).  In  real- 
ity these  terms  were  not  all  measured  in  the  centerline-coor- 
dinate system:  Near  the  wall  Ul/  was  computed  from 


- UV 


as  T 


x*  & 


Since  a hot-wire  responds  essentially  to  resultant 
perpendicular  velocity  and  since  very  near  the  wall  this  is 
virtually  parallel  to  the  wall,  the  (JW  thus  computed  is  actual- 
ly in  the  wall-coordinate  system.  Therefore,  the  second  pro- 
duction term  is  negligible,  and  the  production  is  essentially 

— ui?  Am  . 

The  value  of  i U*  — [/■*)  , ^ 1”  the  centerline  coor- 

i w 


in  the  centerline  coor- 


dinate system  is  shown  in  figure  24. 


"Dissipation  Terms" 
Three  of  the  nine  terms  in 


/ju;  )f  jui  ) 
>dXi<A<jx»/ 


were  measured  using  Taylor's  hypothesis  (ref.  25, 
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following  the  method  of  Townsend  (ref.  26).  The  measured 
terms  are  "T  71  , 7 7 0 , and 

/d_uf  /'JlM2 

(Jx7  gx  J u 

An  estimate  of  the  validity  of  Taylor's  hypothesis  was  made 
following  the  method  of  Ldn  (ref.  -?).  In  a "hear  flow  he  requires 
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Item  (1)  is  valid  with  the  exception  of  very  near  the  wall. 
Item  (2)  would  be  expected  to  hold  for  wave  numbers  larger 
than  k(  (say)  where 

,U  >> 

For  a logarithmic  mean-velocity  distribution 


For  the  present  mean-velocity  distribution  this  reduces  to 

k,b  >>  € s 

where  € has  the  value 
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Figure  19  shows  that  the  greatest  portion  of  > 

which  is  somewhat  comparable  to  the  dissipation  spectrum , is 
at  higher  wave  number  than  Q , except  for  ti  e station  ,/fo  = 
0.01.  In  the  latter  case  ef-io)  bisects  the  area.  With  locally 
isotropic  dissipation  the  dissipation  spectrum  is  actually  the 
second  moment  of  the  three  -dimensional  spectrum,  which,  for 
isotropic  turbulence,  shifts  the  c.  g.  to  higher  wave  numbers. 

In  the  present  experiment  the  values  of  the  "dissipation"  term 
were  most  in  doubt  in  this  region  near  the  wall  due  to  the  large 
magnitude  of  the  terms  which  were  not  measured.  The  estimat- 
ing procedures  are  discussed  below. 
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Having  measured  three  terms  in  W,  the  problem  was 
to  estimate  the  other  six.  Since  Laufer  had  measured  the  same 
three  terms  plus  two  of  the  remaining  terms  in  the  pipe  (ref.  6), 
his  data  were  employed,  on  the  argument  that  very  near  the  wall, 
the  boundary  curvature  should  have  little  effect. 


The  first  step  was  to  select  a suitable  basis  for  estimat- 
ing the  unmeasured  terms  from  those  measured.  The  procedure 
was  to  divide  the  channel  into  two  regions:  (1)0./  a y a I , 
where  b was  used  as  the  characteristic  dimension,  and  (2) 


C £=. 


6 01 


where 


was  selected  as  the  character- 


t i.. 


istic  dimension. 


In  the  center  region  one  would  expect 


X fan') 
I (i  v J 


to  be  a pertinent  parameter.  The  F&  distribution  is  shown  in 
figure  19.  The  values  are  sufficiently  large  that  one  would  expect 
the  dissipation  terms  to  receive  tn^ir  principal  contribution  from 
the  locally  isotropic  fluctuations.  The  values  fall  between 

those  of  the  two  pipe  flows  of  -50,  000  and  500,000.  The 
ratio  between  terms  that  we  did  not  measure  and  those  vve  did 
measure  was  essentially  the  same  in  the  50,  000  and  500.  000 
flows,  so  it  seems  reasonable  to  guess  that  the  same  ratio  would 
apply  here.  The  ratios  obtained  from  the  former  were  used. 


In  the  wall  region  an  assumption  of  similarity  was  used. 
Specifically  we  assumed  that  the  ratio  between  the  unmeasured 
and  measured  terms  was  the  same  function  of  as  in  the 

y 

= 50,  000  pipe.  It  will  be  shown  later  that  some  of  the  other 
turbulent  terms  (e.  g.  ) differ  in  this  region  by  as  much  as 

30  percent  between  the  two  flows.  Since  the  measured  terms  are 
considerably  smaller  than  the  unmeasured  terms  the  local  error 
in  W would  be  large. 


The 

compared. 


next  step  was  to  select  the  form  of  the  terms  to  be 
The  microscales  were  measured  by 
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and  are  compared  in  figure  21.  The  latter  two  agree  very 
well  in  the  center  of  the  channel  , but  they  are  significantly 
smaller  than  , in  particular  near  the  wall.  The  three 

approach  each  other  near  the  center  of  the  channel.  This 
difference  is  expected  since  the  dissipation  terms  (the  deriva- 
tive.?-) essentially  repre  sent^the  high  wave-number  portion, 
wheraa  6<e  , t/*  , and  U T*  receive  their  main  contr  ibution 

from  the  non-isotropic  "feeding"  range. 


Therefore  the  "dissipation"  terms  were  compared 
directly,  figure  22.  The  agreement  is  good.  Therefore  the 
ratios  of  the  dissipation  terms  measured  by  Laufer  were  used 


to  obtain  directly 


and 


(fcf 


The  remaining  four  terms  were  obtained,  following 
reference  6,  by  assuming  that  the  derivatives  in  the  y and  z 
directions  individually  are  equal. 


By  this  procedure  W was  computed  to  be  within  10  percent 
of  the  other  terms  in  the  integral  energy  balance.  It  should 
be  mentioned  that  the  naive  use  of  the  pure  isotropic  relation 


W — 


,,  /d«  t* 
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gave  a value  of  only  60  percent  of  the  remaining  energy  balance. 
However  considerable  question  exists  concerning  the  unmeasured 
terms,  so  an  important  item  in  further  work  would  be  to  measure 


at  least  the  terms 


and 


m 
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DISCUSSION 


The  turbulent  xnotion  will  be  examined  from  three 
principal  points  of  view:  (1)  by  the  distribution  of  energy 
through  the  spectrum,  (2)  through  a comparison  of  the  various 
terms  entering  the  turbulent  energy  balance  and  (3)  in  the 
light  of  the  various  similarity  hypotheses.  Each  treats  the 
motion  from  essentially  a different  aspect  so  they  will  be  dis- 
cussed separately. 


Throughout  the  discussion  emphasis  will  be  placed  on 
the  effect  of  the  adverse  pressure  gradient  on  the  flow 
already  observed  a significant  retardation  of  the  whole  mean- 
velocity  distribution  in  comparison  to  the  channel  case.  Even  more 
striking  however,  was  the  rise  in  the  absolute  level  of  the  turbu- 
lent energy.  In  fact  the  relative  change  in  the  turbulent  energy 
by  the  pressure  gradient  was  larger  than  on  the  mean  motion. 

Also  important  wa3  the  shift  in  the  location  of  the  maximum  shear 
stress  away  from  the  wall  to  ~ 0.  35.  We  will  now 


examine  the  other  effects. 


Spectrum  Analysis 

I Figure  16  reveals  that  a large  portion  (5  < k,b  < 30)  of 
the  u -spectrum  in  the  center  of  the  channel  varys  as  the  -5/3 
power.  Naively  this  could  imply,  following  Kolmogoroff  (ref.  27), 
that  local  isotropy  in  this  flow  extends  down  to  these  wave  numbers 
which  fall  below  the  dissipation  range.  But  the  surprising  point 
is  that  the  - 5/3  power  extends  down  to  wave  lengths  equal  to 
half  the  channel  width,  i.  e.  where  Ztt  = b.  It  is  quite  improb- 

k, 

able  that  the  fluctations  of  this  wave  length  are  isotropic.  A 
similar  situation  was  noted  by  Lin  (ref.  26)  in  the  turbulent 
boundary-layer  measurements  of  Klebanoff  and  Diehl  (ref.  8). 

The  wave  lengths  at  the  low  end  of  the  -5/3  power  region  were 
in  the  order  of  the  bounJary  layer  thickness.  An  equivalent 
situation  exists  in  the  pipe  (ref.  6). 


Searching  for  indications  of  a true  "inertial  subrange"  of 
finite  spectral  extent  we  must  find  that  the  upper  limit  of  the  -5/3 
range  lies  at  a kt  large  compared  with  2tt  . Figure  16  does 

T~ 


show  such  an  inequality. 
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Further,  looking  at  the  v'  and  w'  spectra,  wc  note  that 
they  both  have  -5/3  power  regions  that  do  not  extend  to  as  low 
wave  numbers.  In  tact,  they  are  so  narrow  that  we  cannot 
vouch  lor  their  existence.  l he  same  trend  is  evident  lor  the 
pipe  spectra  (ref.  6).  Evidently  other  factors  have  extended 
the  low  end  of  the  -5/3  power  region  of  the  u-spectrum  (at 
y‘/b  = 0.  1 there  is  a considerable  -5/3  power  region  in  Fu(k, ) 
and  for  Fv(k,)  and  Fw(ki)  it  is  non-existent.  ) However,  the 
u',  v'  and  w!  spectra  do  all  "drop  away"  from  the  -5/3  slope 
at  about  the  same  k,,  a possible  indication  of  local  isotropy 
for  wave  numbers  higher  than  this  one.  The  corresponding 
wave  length,  2 rr  , is  roughly  1 of  the  channel  width.  The 

~ k~  10" 


exister.ee  oi  isotropy  in  the  dissipative 
already  been  pointed  out  by  noting  that 


are  v ery  nearly  in  agr 


spectral  region  has 


cement  with  the  isotropic 


• aw  s . 


One  condition  of  Kolmogoroff' s analysis  to  be  fulfilled 
in  order  to  have  a -5/3  power  region  is  that  in  addition  to  being 
isotropic  the  so-called  inertial  subrange  be  essentially  free  of 
dissipative  effects.  Figure  19  shows  that  75  percent  of  the 
contribution  to 

(°°t 

\ k F (k  ) dk,  at  b = 1.  00 

\ ' u 1 

is  from  wave  numbers  greater  than  these  of  the  -5/3  power 
region.  As  before  it  is  expected  that  even  a larger  portion  of  the 
three-dimensional  spectrum,  v.’hose  second  moment  give3  the 
dissipation  for  isotropic  turbulence  is  above  the  -5/3  power 
region.  Therefore  the  condition  appears  to  be  fulfilled  here. 

The  Kolmogoroff  condition  was  not  examined  at  the  point  nearer 
the  wall  in  view  of  the  negligible  -5/3  power  range  in  their  v1 
ana  Vv  spccci*a. 


The  u'  spectrum  at  the  center  of  the  channel  exhibits  the  -7 
power  range  at  the  high  wave  number  end.  However,  the  other 
spectra  in  the  same  k(b  range  have  slopes  with  a somewhat  more 
positive  Irenu  lOWiii'uS  tii£  wall. 
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In  general  the  shapes  of  the  measured  portion  of  the 
spectra  are  similar  to  the  pipe  spectra  (ref.  6).  The  diffuser 
spectra  do  extend  further  in  the  -7  power  range,  but  the  meas- 
urements were  not  able  to  be  carried  as  far  into  the  low  wave 
number  end  because  the  air  speed  was  much  less.  Inasmuch 
as  the  spectra  in  figure  16  tend  to  be  predominantly  isotropic, 
except  near  the  wall  where  the  viscous  effects  predominate,  any 
influences  of  the  pressure  gradient  or  diverging  walls  would  be 
expected  to  manifest  themselves  in  the  low  wave  number  end 
that  was  not  measured. 


ii*v kdsniul.y  the  cha^act^^  i c w f . . . ^ r. g ... u o . i,.iu  .ui  .iuiCui 

motion  being  so  near  the  overall  dimensions  of  the  whole  turbulent 
shear  region  seems  to  be  inherent  in  the  turbulent  shear  flows. 

For  example,  the  integral  scales  in  the  boundary-layer  are  not 
an  order  of  magnitude  less  than  the  boundary  layer  thichness,  and  the 
velocity  fluctuations  are  found  to  be  correlated  between  points  all 
the  way  across  the  layer  (ref.  9).  The  longitudinal  integral  scale  was 
measured  in  the  present  experiment  by  comparing  the  uncompensated 
and  compensated  responses  of  a hot  wire  by  the  method  described 
in  reference  5.  The  distribution  is  given  in  figure  23.  L^.  is  nearly 

everywhere  equal  to  b,  even  at  y'/b  = 0.  1.  The  values  agree  within 
the  measurement  accuracy  to  Laufer's  values  at  the  same  R,  . 

This  discussion  emphasizes  the  inherent  weakness  of  the  turbulence 
■'theories"  (e.  g.  mixing-length)  which  attempt  to  relate  the 
turbulent  motion  and  the  local  velocity  gradients. 


It  is  well  established  that  the  turbulence  Reynolds  number 
is  a primary  parameter  of  isotropic  turbulence.  It  is  also 
generally  believed  to  be  of  great  importance  in  turbulent  shear 
flow,  especially  in  connection  with  the  relative  extent  of  the  locally 
isotropic  part  of  the  spectrum.  For  example,  for  sufficiently 
high  mx  it  is  expected  that  the  dissipation  will  occur  chiefly  in  the 
isotropic  range;  for  still  higher  , Kolmogoroff's  theory 
anticipates  the  existence  of  an  extended  (-5/3  power)  "inertial 
subrange". 


For  given  boundary  conditions  should  be  uniquely 
related  to  RK,  the  overall  mean  flow  Reynolds  number.  There- 
fore, in  seeking  the  effects  of  changing  boundary  conditions 
(changing  the  channel  half-angle  from  0°  to  1°  ) it  is  important 
to  compare  ' s at  the  same  value  of  R^,  or  to  compare  flows 

with  equal  Ax's  and  different  R^.  In  figure  20  it  is  evident  that 
the  = 61,600  channel  and  the  present  experiment  have 

identical  distributions  over  mo3t  of  the  channel.  The  microscales 
were  computed  in  both  experiments  by 
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which  essentially  compares  the  contributions  from  the  high 
and  low  wave  numbers.  The  Ay  uisii  iuutiow  foi  the  two 
cases  are  in  good  agreement,  figure  21. 


Turbulent  Energy  Balance 

The  conclusions  that  may  be  drawn  from  the  turbulent- 
energy  balance  are  somewhat  restricted  by  the  method  employed 
to  estimate  the  "dissipation11  terms.  Furthermore,  as 
mentioned,  the  true  dissipation  includes  77  \/  i •!  which 

idu<  V <ju«\ 


are  zero  in  the  isotropic  case,  but  not  in  general.  However, 
they  probably  are  small  except  near  the  wall  where  even  the 
value  of  71  i ~ i “ is  very  uncertain.  Nevertheless 
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several  interesting  conclusions  can  be  drawn  from  the  energy 
balance,  shown  in  figures  24  and  25. 


The  flow  field  seems  to  divide  into  four  regions: 

y 1 

1.  The  center  of  the  diffuser,  0.  8 < -j—  ^ 1 where 

a net  amount  of  energy  is  fed  in  by  turbulent  convection  from 

the  low  y'  side  and  by  mean  flow  convection  from  the  upstream  side. 

This  latter  effect  coneitiutes  a difference  from  the  case  of  the 
parallel  wall  channel. 

2, qi  ^ y 1 ^ q g This  region  differs  from  the 

b 

channel  case  in  that  here  (a)  the  production  of  energy  is  mostly 
balanced  locally  by  the  dissipation  and  (b)  the  energy  convected 
into  the  region  by  mean  flow  is  very  nearly  offset  by  the 
turbulent  diffusion  and  pressure  work  "outward". 

i 
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.3.  0.  05<  - ^ 0.1.  In  this  region  similar  to  the  pipe 

b 

(ref.  6),  the  production  and  dissipation  are  in  balance  and  the 
remaining  terms  are  relatively  smaller.  This  means  that  the 
energy  produced  here  is  locally  dissipated,  a result  consistent 
with  the  render  assumptions  of  the  mixing-length  theories.  Hence 
this  may  be  regarded  as  consistent  with  the  occurrence  of  a 
logarithmic  mean- velocity  distribution  in  this  region,  to  be  shown 
later. 


& 


4.  Wall  region.  As  in  the  pipe 
a maximum  at  about  n.h  m - n. 


J M 


the  production  reaches 
just  outside  the  laminar 


sublayer.  If  we  assume  a qualitative  resemblance  to  the  wall 
region  of  the  pipe,  it  can  be  inferred  that  the  energy  is  transported 
away  by  the  turbulent  diffusion  and,  to  a lesser  extent,  b y the 
"viscous  work"  S'  the  true  viscous  work  includes  /Jlii  \/  ^ 

^ (j  XtJlaWj  j 

However,  in  contrast  to  the  pipe  case,  the  dissipation  in  this  region 
exceeds  the  production,  so  we  must  conclude  that  the  feeding  of 
energy  into  this  region  by  pressure  work  (Laufer's  "pressure 
diffusion")  exceeds  the  two  diffusion  effects  which  tend  to  drain 
off  the  energy.  This  difference  between  the  wall  region  of  a pipe 
and  diffuser  may  be  important  if  correct,  since  it  is  widely 
assumed  that  this  small  region  near  the  wall  is  "unaware"  of  the 
flow  in  the  large,  that  is  the  overall  Reynolds  number  and  pressure 
gradient.  Therefore,  a more  detailed  measurement  of  the 
dissipation  in  this  should  be  a part  of  further  work  in  order  to 
improve  the  energy  balance  and  thus  learn  whether  this  effect  of 
an  adverse  pressure  gradient  does  exist. 


Dimensional  Reasoning 

It  is  well-known  that  the  turbulent  boundary  layer  can 
be  separated  for  analytical  purposes  into  two  regions:  one  compris- 
ing all  the  flow  field  except  a thin  layer  nearest  the  wall,  which 
has  as  a characteristic  length  the  boundary-layer  thickness,  and 
velocity  xx^.  ; and  the  other  being  the  wail  region  with  a character- 
istic length  V , and  velocity  u^.  . These  correlations  will 

W* 

be  examined  for  the  channel  and  diffuser. 

1,  Center  of  the  channel.  The  mean  velocity  is  compared 
on  the  basis  of  the  velocity-defect  law  in  figure  26.  For  the  turbu- 
lent b ou  d u r y layer  this  representation  of  the  mean- velocity 
distribution  has  been  empirically  determined  to  be  independent  of 
the  Reynolds  number  and  surface  roughness,  but  a function  of  the 
pressure  gradient.  Since  all  the  diffuser  tests  were  conducted  at 
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a constant  Reynolds  number  a oingle  curve  is  expected.  The 
trend  of  the  channel  and  diffuser  data  with  an  adverse  pressure 
gradient  is  similar  to  the  boundary  layer.  Of  course,  the 
boundary  layer  and  diffuser  data  are  expected  to  differ  from  each 
other  in  this  region  and  only  the  trends  are  compared.  The 
pressure  - gradient  parameter,  ^ dP  wa  a 7.  7 for  the 


diffuser,  and 


jL  if 


7“ 


was  12.  2-14.  8 for  the  boundary  layer. 


The  turbulence  level  was  compared  on  a similar  basis, 
figure  27.  This  emphasizes  the  larger  effect  of  the  pressure 
gradient  on  tire  turbulent  fluctuations  than  on  the  mean  flow. 


2.  Near  the  wall  the  results  are  considerably  more 
interesting.  The  mean  velocity  is  compared  on  the  basis  of 
Prandtl' s wall- or oximity  law  (ref.  29;  in  figure  28  using  semi- 
logarithmic  coordinates  to  demonstrate  the  log-law  region. 

4 The  log  law  extends  over  a wide  region  up  to  y^t  = 150. 

The  difference  in  displacement  of  the  various  profiles  m this  region 
is  probably  not  significant,  in  fact  all  of  the  points  could  be 
, collapsed  to  ^-re  curve  by  admitting  a 10  percent  error  ir  u*. 

Therefore  the  mean  velocity  distribution  could  be  said  to  fulfill 
tne  similarity  hypothesis  out  to  y*  150. 
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The  turbulence  level  story,  however,  is  quite  different, 
figure  29.  Laufe  r found  good  agreement  between  the  channel 
(ref.  5)  and  the  pipe  (ref.  6),  with  some  effect  of  Reynolds 
number  for  y+  > 30.  But  the  effect  of  the  pressure  gradient 

is  evidently  much  larger,  and  it  is  not  attributable  to  instrument 
errors.  In  fact  a correction  in  u*  to  bring  the  channel  and 
diffuser  data  into  better  agreement  in  U , figure  28.  actually 
increases  the  discrepancy  in  u',  figure  29.  Some  unpublished 
data  was  obtained  from  Dr.  Clauser's  group*  which  had  measured 
this  term  in  the  turbulent  boundary  layer  in  an  adverse  pressure 
gradient.  The  values  at  four  stations  were  all  grouped  about  the 
«'  = 1 deg.  curve.  The  pressure  gradient  based  on  the 
characteristic  quantities  in  this  region  was  somewhat  smaller. 

The  values  of  V dP  were: 

piP  "3*  , 

Turb.  b.l.  5.4  x 10"^  to  2.  9 x 10‘3 

Diffuser  6.8  x 10 


^Aeronautics  Department,  Johns  Hopkins  University. 


An  examination  was  made  to  ascertain  wnether  u*  is 
a characteristic  velocity  of  the  turbulent  region  since  'V 
changes  across  the  laminar  sublayer.  But  the  variation  in 

L Was  omy  D percent.  lllouixiC  iciik  . O aCCvuii.  the  d.  — 

The  question  of  local  similarity  was  examined  further  by  comparing 
U,‘_r  , which  eliminates  the  characteristic  velocity.  This 

air7 

is  shown  in  figure  30.  The  dashed  line  represents  the  diffuser 
values  in  the  wall-coordinate  systems.  Of  course,  u#  remains 
in  the  abcissa  but  the  difference  between  the  curves  is  so  great 
that  no  reasonable  shift  by  adopting  a new  characteristic 
velocity  could  collapse  the  two  to  one. 

Another  approach  was  then  taken.  Since  only  the 
longitudinal  motion  faces  the  mean  adverse-pressure  gradient 
the  pressure  effects  may  have  a preferred  direction.  The  ratio 
of  u'  to  v'  was  compared  between  the  pipe  flow  and  the  diffuser 
flow,  figure  30,  and  the  agreement  in  the  wall- proximity  region 
is  surprising.  The  conclusion  to  he  drawn  evidently  is  that  the 
pressure  gradient  docs  not  affect  the  turbulence  levels  near  the 
wall,  but  only  their  correlation. 

The  distributions  of  u'  at  equal  Reynolds  numbers  are 

v ' 

also  compared  in  the  center  of  the  channel,  pipe,  and  diffuser  in  fig. 
30.  In  this  region  the  difference  is  significant.  Evidently  tne 
pressure  gradient  does  manifest  itself  by  favoring  the  fluctuations 
in  the  direction  of  the  gradient,  perhaps  in  an  indirect  way  through 
the  energy  balance.  For  example,  the  production  feeds  almost 
purely  into  u-fluctuations , while  turbulent  diffusion  acts  on  all 
three  components. 

Several  oscilloscope  traces  were  made  of  u(t)  within  and 
slightly  outside  the  laminar  sublayer,  figure  31.  The  fluctuations 
inside  the  sublayer  are  interesting.  They  appear  to  display  an 
intermittancy  suggestive  of  the  outer  regions  of  a boundary  layer. 
An  important  difference  is  the  absence  of  the  sharp  demarcations 
between  the  turbulent,  and  in  this  case,  laminar  regions,  that 
is  evident  in  boundary-layer  intermittancy.  Another  difference 
is  the  fact  that  the  turbulent  regions  near  the  wall  are  generally 
traveling  with  excess  velocities,  whereas  in  the  boundary  layer 
they  tend  to  move  slower. 
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Further  Investigations 

It  seerns  apparent  following  the  present  investigation 
that  further  work  should  pursue  the  following  course. 

1.  Measure  several  of  the  dissipation  terms  in  the 
direction  normal  to  the  main  flow,  in  particular  near  the  region 
of  maximum  production. 

2.  Measure  the  triple  -correlation  distribution  to 
identify  the  pressure-  and  velocity-diffusion  terms.  Laufer'a 

data  indicate  that  both  of  these  effects  are  very  large  near  the 
maximum-production  point. 

3.  Look  for  increased  skewness  in  some  turbulence 
properties  - a possibility  suggested  by  the  unsymmetrizing  effect 
of  a large  static-pressure  gradient. 

4.  Increase  the  angle  of  divergence  to  find  the  effect  of 
the  pressure  gradient  as  the  flow  approaches  separation.  In 
particular  to  study  the  similarity  considerations  in  the  wall 
proximity. 
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APPENDIX  I 


ESTABLISHMENT  OF  THE  EXPERIMENTAL  FLOW 
Setting  Up  The  Test  Conditions. 


The  most  difficult  phase  of  setting  up  the  test  con- 
ditions was  making  the  mean  flow  reasonably  two-dimen- 
sional. This  is  apparently  a problem  characteristic  of 
adverse-pressure  gradient  flows  although  it  has  only  recently 
received  attention  in  the  literature.  Ludwieg  and  Tiiimann 
(ref.  10  ) found  that  by  including  the  measured  secondary  flow 
in  a boundary  layer  in  an  adverse-pressure  gradient  the  value 
of  Cf  computed  by  von  Karman's  momentum  equation  was 
reduced  by  40  per  cent.  Clauser  (ref.ll  ) found  even  larger 
discrepancies  with  very  strong  adverse  pressure  gradients. 

He  indicated  the  importance  of  and  how  it  far  outweighed 

the  terms  mentioned  by  others. 

In  the  present  experiment  the  steps  taken  to  assure 
longitudinal-similarity  made  it  difficult  to  attain  two-dimension- 
ality. The  parallel- wall  approach  channel  was  added  ahead  of 
the  test  section  to  give  the  flow  a start  towards  reaching  the 
fully- developed  state  and  to  give  a previously  studied  flow  field 
in  which  measurements  could  be  checked.  But  inevitably  non- 
uniformities develop  in  the  vertical  direction  by  the  time  the 
stream  has  passed  through  the  parallel  channel,  and  these  are 
amplified  in  the  adverse  pressure  gradient.  The  same  phenomenon 
occurs  when  a boundary  layer  grows  on  a flat  plate  and  enters  an 
adverse  pressure  gradient. 

Some  physical  feeling  for  the  phenomenon  can  be  reached 
by  the  following  rather  rough  analysis  of  the  flow  in  the  present 
diffuser.  Consider  the  stream  at  two  elevations:  the  midplane 
( 2*0  ) and  2 1/2  in  from  the  ceiling  (2*  3! in  ).  The  ratio  of  the 
maximum  dynamic  pressure  in  the  upper  profile  to  the  maximum 
dynamic  pressure  in  the  midplane,  both  at  the  same  0C  station, 
will  be  called  the  dynamic-pressure  ratio,  Pf  . fr  is  unity  if 
the  flow  is  two-dimensional.  The  measured  values  of  Pr  are 
tabulated  in  figure  32. 


1 


46 


r 


\ 


i 


Let  us  naively  invoke  the  Bernoulli  equation  for  the 
flow.  We  will  consider  the  maximum  dynamic  pressure 
{presumably  at  ‘j=C  ) for  thi«  approximation.  When  the 
walls  of  the  parallel  section  were  smooth  the  value  of  /"£. 

6 in.  ahead  of  the  diffuser  inlet  was  0.  96:  a velocity  defect 
of  only  2 per  cent.  With  ?.<t  - 2 degrees  the  width  at  station 
2 is  2.  9 times  the  width  in  the  approach  channel,  so  a 
Bernoulli  expansion  at  the  midplune  would  convert  88  per 
cent  of  the  original  dynamic  pressure  into  static  pressure, 

12  per  cent  remaining  in  the  dynamic  pressure.  The  upper 
plane  must  convert  the  same  amount  of  energy  into  static 

pressure  because  -2.C.  is  small.  But  it  only  had  0.  96  as 

' d £ 


much  to  start  with, 


so  Bernoulli-wise 


Pr 


096-088 

would  be _ 

1.00  — u8B 


or  0.  67,  the  value  measured.  The  33  per  cent  defect  re- 
sult ed  from  the  small  4 per  cent  defect  at  the  inlet.  At 
station  5 the  dynamic  pressure  at  the  top  would  be  negative 
by  Bernoulli's  equation,  so  momentum  had  evidently  been 
transferred  away  from  the  middle  region,  a condition  to  be 
avoided.  (Separation  did  not  take  place,  in  fact,  anywhere 
in  the  diffuser  at  2 «C  = 2 degrees.  It  did,  however,  at  3 
degrees) 


Nikuradse's  (ref.  13  ) furthest  downstream  station 

for  2<C=  2 degrees  wae^=  34.  7,  roughly  half  that  of  station 
2,  and  he  measured//.  = 0.985.  In  his  experiment  the  flow 
entered  the  diffuser  through  a contraction  so  the  velocity 
distribution  was  initially  uniform.  The  question  arises  as 
to  whether  he  had  a fully  shearing  flow,  or  just  boundary 
layers. 

The  measurements  of  Pf.  in  the  approach  channel 
shewed  that  the  two-dimensionality  troubles  were  located 
in  the  upstream  portion.  Apparently  the  following  was 
happening:  at  the  contraction  from  the  plenum  chamber  the 
floor  and  ceiling  are  plane,  the  contraction  being  all  in  the 
Bide  walls,  therefore  the  boundary  layers  on  the  floor  and 


■"■at* 
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ceiling  tend  to  be  thicker.  Furthermore,  the  interaction 
of  the  corners  retards  the  boundary  layers  on  the  ends  even 
more.  Therefore  transition  started  first  at  the  floor  and 
ceiling  and  further  downstream  at  the  midplane.  Since  the 
coefficient  of  friction  of  a boundary  layer  roughly  doubles 
after  passing  through  transition,  the  momentum  was  removed 
at  a faster  rate  near  the  floor  and  ceiling  until  the  flow  at  the 
midplane  passed  through  transition. 

The  transition  zone  was  measured  at  several  eleva- 
tions by  dragging  an  impact  probe  along  the  surface.  The 
method  is  similar  to  that  employed  by  Klebanoff  and  Diehl 
(ref.  8 ) with  the  exception  that  transition  was  determined  by 
the  total  pressure  in  the  present  experiment  instead  of  the  im- 
pact pressure.  At  the  midplane  transition  started  about  12in. 
downstream  from  the  point  where  the  inlet  contraction  ends 
and  the  parallel  wall  begins,  and  had  not  ended  by  20in.  from 
the  inlet.  At  2 l/2in.  from  the  ceiling  transition  took  place- 
within  4in.  downstream  of  the  inlet. 

Two  methods  were  tried  to  improve  the  two-dimension- 
ality of  the  flow  leaving  the  approach  channel.  In  the  first  plan, 
which  was  the  more  conservative,  a row  of  tacks  was  imbedded 
in  each  wall  12in.  from  the  inlet.  The  heads  were  0.  22in,  dia. 
znii  0,  02 Sir.  thick  and  projected  entirely  above  the  surface  the 
intent  being  to  remove  momentum  by  the  profile  drag.  The 
interval  between  tacks  was  increased  near  the  floor  and  ceiling 
to  reduce  the  momentum  removal.  It  turned  out  that  the  improve- 
ment of  the  two-dimensionality  was  negligible. 

Thereupon  we  decided  to  utilize  the  difference  between 
the  laminar  and  turbulent  Cf  by  artificially  causing  transition 
further  upstream  along  the  middle  region.  At  this  time  the 
problem  was  discussed  with  Mr.  Phillip  Eisenberg  of  the  David 
Taylor  Model  Basin.  * In  particular  he  indicated  that  single 


* 

Currently  at  the  Office  of  Naval  Research  (Head,  Mechanics 
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roughness  elements  had  been  used  succes siuily  to  promote 
transition  on  ship  models,  anci  he  directed  us  to  the  litera- 
ture on  this  work  and  referred  us  to  the  research  that  has 
been  done  by  the  David  Taylor  Model  Basin,  One  advantage 
of  single  roughness  elements,  as  he  painted  out,  lo  their 
small  profile  drag.  In  the  channel  they  have  an  additional 
merit.  The  stream  contracts  as  it  passes  over  the  tripper 
and  the  pressure  gradient  tends  to  become  more  favorable 
since  Bernoulli’s  equation  applies  to  the  stream  core.  If 
the  transition  device  is  continuous;  such  as  a wire,  the  favor- 
able pressure  gradient  could  be  so  large  that  the  stream  could 
conceivably  reattack  and  remain  laminar. 

We  discussed  the  problem  with  Mr.  Marshall  Tulin  of 
the  David  Taylor  Model  Basin,  and  his  experience  with  short 
cylinders  normal  to  the  surface  was  very  helpful  to  us.  Also 
the  influence  of  Reynolds  number  {ref.  30  ) for  other  trippers 
could  be  interpreted  to  our  case.  Therefore  we  decided  as 
the  second  plan  to  place  a row  of  tacka  on  each  wall  at  a distance 
of  3 1/4  in.  from  the  inlet  at  the  midplane  and  curving  further 
downstream  near  the  floor  and  ceiling.  The  initial  results  were 
surprisingly  good.  Following  several  hours  running  we  dis- 
covered, quite  by  accident,  that  the  dust  formed  on  the  surface  in 
such  a manner  that  when  a light  was  placed  nearly  flat  to  the  sur- 
face the  laminar  regions  were  dark  and  the  turbulent  regions  light.  * 
The  dust  patterns  revealed  that  some  of  the  tacks  had  not  caused 
transition.  We  found  these  tacks  did  not  protrude  as  far  from  the 
surface.  They  were  pulled  further  out  whereupon  all  the  tacks  had 
turbulent  boundary  layers  in  their  wakes. 


* 

No  explanation  has  been  given,  in  fact  one  would  expect 
the  opposite  to  take  place.  When  the  surface  was  wiped  it  ail 
became  dark,  so  it  is  clear  that  dust  had  settled. 
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We  then  remeasured  at  the  cutlet  of  the  approach 
section;  at  the  top  and  bottom  it  improved  to  1.  00  and  0.  985 
respectively.  The  two -dimensionality  in  the  test  section  was 
also  significantly  better  and  the  very  low  frequency  fluctuations 
that  interfered  with  velocity  measurements  were  considerably 
reduced.  A ?mal)  Reynolds  number  effect  and  a difference  when 
the  number  of  fans  were  changed  was  observed.  This  was  not 
traced  further.  Some  effort  was  spent  to  raise  the  dynamic  pres- 
sure at  the  floor  but  with  little  success.  These  values  were  there- 
fore accepted.  The  characteristics  of  the  downstream  flow  were 
measured  in  some  detail,  so  they  will  now  be  examined. 

Evaluation  Of  The  Test  Conditions 

The  purpose  of  this  section  is  to  determine  how  well  the 
wind  tunnel  satisfied  the  assumptions  made  in  the  Preliminary  Analy- 
sis. We  will  show  how  the  integral  equations  derived  in  the  Pre- 
liminary Analysis  section,  can  be  used  to  make  definite  numerical 
statements  about  the  individual  conditions. 


In  principle  the  plan  will  be  to  examine  each  term  that  was 
eliminated  in  the  analysis  by  virtue  of  the  assumptions  and  to  deter- 
mine whether  it  actually  is  as  small  as  the  analysis  requires.  This 
will  eliminate  the  possibility  of  two  cancelling  affects. 


The  examination  of  the  basic  assumption  concerning  contin- 
uity, the  Navier-Stokes  equations,  and  the  Reynolds  postulate  uo  ob- 
tain equations  (1),  (2),  and  (3)  is  more  sophisticated  than  we  are 
attempting  here  so  we  will  proceed  from  these  equations.  In  expan- 
ded form  equations  (1)  and  (2)  are 


aO  + av  . $w s 0 
3x  ay  a* 


d') 

(2 -a’  ) 


u ax  ay  az  f d*  l<)z  1 

> jf  * 5 f ♦ * f,-  - lf  af  ‘ af ' ' Wj 


r-^m-  -«,|P  py«S*’ 
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A complete  discussion  of  turbulence  assumption  (d) 
and  (e),  and,  more  particularly,  hov/  they  were  used  is  quite 
involved  and  would  divert  attention  from  the  more  important 
assumptions.  Also  the  main  items  were  covered  in  the  sec- 
tion on  Preliminary  Analysis.  The  question  herq_js  whether 

the  flow  satisfied  the  conditions:  (a)  W-0  , (b)  au.o 

02 

for  all  terms,  the  boundary-layer  approximation,  and  longi- 
tudinal similarity. 

From  the  boundary-layer  approximation  we  deduced 
that  \I«U  • Measurements  at  tne  ?.n  1 d p ian e oi  station  5 \vith 
the  x-meter  showed  this  to  be  the  case.  Furthermore  was 
measured  by  the  x-meter  and  claw  probe  at  the  same  location 
and  it  was  significantly  smaller  than  t/_  , figure  33  But  the 

derivatives  , and  5v  &{/  must  be  carried  since 

d*_  _ 62 

is  small.  Since  W « W , then  each  term  in  equation 
{2-c  -)  is  smaller  than  its  counterpart  in  equation  (2-a*).  But 
it  is  possible  that  in  the  aggregate  they  may  interact  through 
the  pressure  upon  equation  (2-a*).  It  will  be  demonstrated  later 
by  measurements  that  d]P  was  sufficiently  small,  sc  we  will 

az  . i 

use  this  fact.  Pressure  interaction  from  equation  (Z-b  ) is  not 

likely,  except  through  Vt  as  mentioned  in  equation  (5b);  because 

V was  less  than  one-fourth  of  W . Therefore  equations 

(1*).  (2-a*),  (2-b*),  and  (2-c*)  reduce  to 


a 0 
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Equation  (!')  can  be  substituted  into  equation  (t  ) giving 

duz  , d\u  v)  . dww)  _ _ l 4_P\  i)  £jj  _ d_My 
dX  dy  * £ P dX  dys  dy 

This  equation  will  now  be  put  into  the  integral  form  to  study 
the  remaining  assumptions.  The  second  term  drops  out 
leaving 

(rf'ju  + Urn<iu  = -t*  dP‘h  i>[0\  -0)  1+fe,) 

jix^Jaz  y P dr  Twy-L  UoKd  ‘ V' 


The  flow  symmetry  was  ch 


ecked  and  I %? ) « f£P) 

l^y/y.o  Jy  /y*fc 


station  5 in  the  midplane.  Figure  13  shows  that  (aav)  was 

y*o 

negligible,  so  the  terms  on  the  right-hand  side  of  the  equation 
reduce  to  the  pressure  term  plus  ^-fiu  . The' critical 
assumptions  are  wrapped  up  in  the  left-hand  side.  We  will 
now  expand  these  terms.  The  first  term  is 


JfP* 


’ - "*  :# 


and  since  Ufa  =0  it  can  be  further  expanded  to 

= u‘fb  & + aM ' _ z,  * 7 

dx  do  dx  J 

X 

w>  - . __  r J_  a Si 

or  non- dimens ionally,  where  p (XI  ZZ  J j 1 

L d?  7?t  d?  -J 
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The  second  term  on  (.he  left  hand  side  is 


J Z* 


U.’[ + z f JOLfit  J’j.l  . fs2-*dv 7 

U d(f)  < J u fa&)  < atm)  y u j 


so  the  complete  integral  equation  in 


If?  + Sh  +2[ V7~%  jf  ff&dV= 

d?  Kb  <3?  J d(f)  1 J u d[f)  - V'J(f)jf  D I 


L¥*  + Iz* 


/°4>*  c/X 

6 + 


Zo 

/>vt 


(IV) 


which  differs  from  equa.tion  (13)  in  the  Preliminary  Analysis 
by  the  terms  on  the  left-hand  side.  These  terms  drop  out 
if  the  assumptions  are  realized.  The  first  term,  1 + , vanishes 
if  the  mean-velocity  profiles  satisfy  longitudinal  similarity.  If 
the  profiles  are  similar  but  fluid  is  fed  into  the  plane  then  the 
second  term,  2 + , is  positive.  Terms  3^,  4 + and  5+  are  the  two- 
dimensionality  effects.  Term  3+  is  the  removal  of  momentum 
due  to  a spreading  of  the  stream  surfaces  above  and  below  the 
plane;  term  4+  incorporates  the  two-dimensionality  of  the  mean- 
velocity  t*  i i n 1 1 3. r it y term  5 ^ £ic  counts  for  h.  vertical  in  ths 

centerline  dynamic  pressure.  At  the  midplane  4+  and  5*  would 
be  expected  to  be  zero  by  vertical  symmetry.  Terms  4+  and  5 + 
are  sensitive  to  the  deviation  of  the  stream  surface  lrom  plane - 


ness. 
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The  separation  of  *he  longitudinal  changes  into  the 
form  of  1+  and  2+  and  the  vertical  into  3',  4+;  and  5+  was 
done  to  facilitate  the  experimental  check.  Analytically  it 
might  be  preferable  to  use  the  mean  velocity  across  a section 
instead  of  ; then  1+  would  check  longitudinal  similarity, 

and  2 + , 3 + , 4+,  and  5"*  would  check  two-dimensionality. 


r..i  f. ii  .v.  „ 

4UA  kAA* 


a"  SUHipti™  ® O ^ analvai  q 


w 


would  like  to  have  1+, , 5^  individually  smaller,  by  say  a 

fafctor  of  10,  than  any  of  the  terms  that  remain,  i.  e.  6+.  ?+.  or  8+. 
At  fit  = 1 degree  the  pressure  gradient,  6+,  and  momentum  gra- 
dient. 7+,  are  nearly  in  balance  and  each  is  about  10  times  larger 
than  8+,  which  is  13  x 10“^.  Therefore  the  requirement  is  that 

1+, , 5+  <<  13  x 10"4.  In  addition  we  should  like  to  avoid  local 

irregularities  in  the  flew.  The  limit  of  integration  in  1 + , 3 + , 4 + , 
and  5^  is  unity  so  we  are  on  the  conservative  side  if  we  ask  that 
the  integrand  everywhere  remain  in  the  same  bounds. 


Term  1+  between  stations  2 and  3 was  l.lxlO^so  it  is 
sufficiently  small,  thus  satisfying  longitudinal  similarity.  Let  us 
examine  the  integrand  in 
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Between  stations  2 and  3 the  requirement  would  be  that 
M < 35  (10“4).  Figure  8 shows  this  is  about  the  du- 
plication of  the  profiles, 

H **■*•».  in  where  a question  arises  in  regard  to  Nikuradses' 
experiment  (reference  13  ).  His  furthest  downstream  station 

was  at  JS  = 34.  7 for  ®C  = 1 degree.  Results  of  pipe  flow  indicate 
that  ^ values  of  50  to  100  are  required  for  uniform  flow  at  the  in- 
let to  become  fully  developed.  Nikuradse  made  additional  velocity 
measurements  at  = 23.  2 to  determine  that  the  flow  was  fully 
developed.  The  above  test  would  require  fa 17  >:  10-4.  or  since  f 
has  a value  of  about  unity  the  velocity  profiles  must  duplicate  to 
0.2%.  This  is  better'  than  customary  accuracy.  Whether  the  pro- 
files were  this  close  is  not  known  since  Nikuradse  did  not  report 
the  up  str  eam  data. 


In  the  present  experiment  term  1 between  stations  2 and  5 
was  - 1.  7 x 1 0~4.  Again  this  is  sufficiently  small.  For  the.  in- 
tegrand to  be  everywhere  small  we  require  Uf  0.  004.  Several 
points  were  larger,  but  the  discrepancy  did  not  appear  serious.  To 
measure  1 + the  precision  manometer  was  needed  as  it  corresponds 
over  a large  region  of  the  channel  to(£z)  0*  004  between  the 

^ tfo/7s«wwt. 


stations.  The  dynamic  pressure  at  station  5 for  Up  is  0.  06  in. 
of  water,  so  the  dynamic  pressure  should  be  resolved  to  5 x 
10”4  in.  ^of  water.  The  final  manometer  v/as  more  accurate 
than  10”**  in.  of  water.,  and  is  discussed  in  Appendix  II.  The 
accuracy  of  measuring  T was  tested  by  periodically  measuring 
U0  during  each  traverse  and  comparing  it  to  the  reference 
pressure.  It  had  a probable  error  (0.675  x standard  deviation) 
of  0.  7 x 10-4  in.  of  water  which  is  satisfactory. 


The  mean-velocity  profile  at  stations  4 and  5 was  slightly 
more  peaked  than  at  the  upstream  stations.  This  can  be  related  to 
term  2.  After  the  peaking  was  observed  & was  measured  and 

. P dx 

between  stations  3 and  5 increased  4 percent.  But  during  each 
velocity  traverse  the  wind-tunnel  speed  had  been  set  so  /ffe  =31,400 
at  the  particular  station.  Therefore  when  the  downstream  traverses' 
were  made  the  upstream  flow  was  at  a lower  , and  thus  more 

peaked,  than  when  the  upstream  traverse  had  been  made.  Laufer's 
data  (ref.  5 ) would  indicate  that  this  would  be  a rise  in  f ath  = 0.  2 


1 
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of  about  0.  005.  This  would  account  for  roughly  one-half  of 
the  change.  This  4 percent  rise  in  Reynolds  number  was 
confirmed  by  a corresponding  increase  in  the  dynamic  pressure 
at  the  inlet  to  the  parallel-wall  approach  section  between  the 
two  runs.  The*  latter  method  would  indicate  an  increase 

between  stations  2 and  3 of  3.  5 percent.  The  net  result  is  that 
2+  between  stations  2 and  5 was  21  x 10"^,  or  roughly  equal  to 
the  turbulence  term,  8+.  It  was  only  slightly  balanced  by  1+. 


Term  3 is  believed  to  be  the  offender.  Unfortunately  it 
is  very  difficult  to  measure  directly:  for  it  to  be  smaller  than 

/q  at  station  5 require t.  that  the  iicw  inclination  differ 

by  less  than  0.  04  degrees  between  the  midplane  and  the  traverse 
points  16  in.  above  or  below  the  midplane.  However,  an  estimate 
of  its  magnitude  can  be  made  from  continuity  of  the  flow. 
Integration  of  equation  (1)  across  the  channel  using  only  the  wall 
no- slip  condition  and  V = O at  ij  = 0 gives 


T+  4-  + 

Terms  4 and  5TT  can  be  neglected,  as  before,  by  vertical 
symmetry.  Between  stations 

l A ^ 1*  tl  A 


10 

10 


and  term  2 


4 

negative  as  expected 


was  11.  9 x 10 
If 


and 

-4 


5 term  1+*  was  -0.  7 x 
Term  3 + + is  then  -11.  2 x 
is  independent  of  7^ 


UJL 

2 


m 


it  would  have  the  value  of  ~13  x 10 


Then  putting  this  back 

in  equation  (13  ) would  indicate  that  term  3+  is  -10  x 10' *. 
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The  following  table  is  a summary  of  the  terms  in  equa- 
tion (13*  ) 


Term 

i+ 

* 

1 

4 

7 x 10 

H 

2 + 

* 

21 

it 

II 

3 + 

- 10 

ti 

..  -i- 

1 1 

b • 

rr 

negligible 

Total 

9 

■ x 10^ 

Term 

6+ 

-100 

X 

O 

1 

II 

7+ 

— 

127 

ri 

II 

Q + 

u 

5 

- 13 

n 

Total 

14 

it 

Nuic;  term  4 was  not  measui  eel,  but  i □ believed  to  be  small. 


It  is  evident  that  terms  ?.+  and  3+ roughly  account  for  the  dis- 
crepancy of  13  x 10"^  in  equation  (13).  The  von  Karman  momen- 
tum equation  for  the  turbulent  boundary  layer  assumes  terms 
equivalent  to  3 , 4,  and  5 are  zero.  In  view  of  the  present 

remarks  it  is  not  surprising  that  in  the  past  the  use  of  the 
von  Karman  equation  to  compute  Xi.  has  led  to  apparent  in- 

rfo 

creases  in  the  wall  shear  stress  during  a retardation,  whereas 
a decrease  is  expected. 


The  net  result  is  that  longitudinal  similarity  existed, 
within  the  requirement  set,  but  two-dimensionality  of  the  mean  flow 
deviated  by  an  amount  equal  to  the  turbulence  term  in  equation  (13  ) 

. The  hope  is  that  since  this  effect  acts  through  the  mean 

/>  On- 


flow it  would  not  appreciably  disturb  the  turbulence.  What  it  may 
amount  to  is  a 10  percent  reduction  in  the  "effective"  **  . One  is 
restrained,  however,  from  concluding  at  this  point  that  the  assump- 
tion of  longitudinal  similarity  for  a fully  two-dimensional  channel 
is  confirmed. 


At  this  point  consideration  was  given  to  diverging  the  floor 
and  ceiling.  However,  this  involves  mechanical  difficulties,  and  in 
view  of  the  sensitivity  of  the  flow  to  very  small  changes  in  this 
divergence  angle,  as  indicated  by  term  3^in  equation  ( 13')  , the 
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decision  was  made  to  measure,^  at  two  stations.  If  this 
were  similar,  as  it  turned  out  to  be,  then  it  would  be  con- 
sidered sufficient  evidence  of  flow  similarity  to  proceed 

• <1  r il. A 
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The  M,'  profile  at  3taticn  2 is  compared  to  the  pro- 
file at  station  5 in  figure  10  . It  compares  very  well.  But 
again  the  statement  can  be  made  more  precise,  this  time  by 
means  of  the  turbulent  energy  equation.  A complete  analysis 
including  all  of  the  assumptions,  as  was  done  for  the  momen- 
tum equation,  will  not  be  undertaken.  But  the  similarity 
assumption  will  be  carried  out.  This  was  done  by  rederiving 
equation  (14)  without  making  the  assumption  of  similarity.  The 
following  additional  term  is  obtained,  d I#  . It  can  be 

2 d% 

expanded  to 


/ 


For  each  of  these  terms  to  be  negligible  it  must  be  small  com- 
pared to  the  smaller  of  the  terms  that  remain,  which  is 


The  latter  term  will  be  estimated 


by  the  isotropic  relation 

c - 

9\ 

(dzJ(3*J 

J 

A 

where  the  values  of 

A 

and 

tZ 

■Z—  at  station  5 are  used. 

Thus 

b 
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only,  then77 
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are  function  of 
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to  be  less  than 


we  require  ~r  and  to  be  less  than 

2 ^ ; d?  *a  . .5 

10  . Between  stations  2 and  5 this  means:  "IT*  *7”  ^ 

'I’ll  lo  time  ina  I I n ■)  o t -!  a ® 


This  was  well  satisfied. 


There  is  one  point  to  be  cleared  up.  The  state- 
ment was  made  that  the  interaction  between  equation  (2  a) 
and  (2  c*)  through  the  pressure  is  small.  This  will  now  be 
demonstrated.  If  equation  ( 2-c ')  is  put  in  the  integral  form, 
similar  to  equation  (lV  ),  the  pressure  term  is 
/ 

f b d P J~> 


JL  *£dr? 
J f>US  6%  * 


The  wall  values  of  AZ:  were  measured  between  the  mid- 

plane and  16  in.  above  and  below  the  midplane  to  estimate  this 
term.  It  is  compared  to  term  6 of  equation  (13  ) in  figure  33 
and  it  stays  well  below  one -tenth  of  the  latter  term.  To  be 
assured  that  dP  does  not  affect  the  turbulence  one  m ight  like 

dt  i r 

this  ;~.teeral  to  be  even  less  than  — — — „ . This  is  not  quite 

10  pus 

satisfied  downstream  of  about  station  4 , but  the  discrepancy 
is  not  believed  to  be  serious. 
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APPENDIX  II 

PRECISION  INCLINED  MANOMETER 
Introduction 

In  order  to  measure  mean  velocities  with  sufficient 
reliability  to  establish  the  existence  of  simple  similarity.it 
is  desirable  to  measure  pressure  with  an  accuracy  of  0.  0005  in. 
of  water  (this  is  justified  in  Appendix  I.  ).  The  same  is  re- 
quired to  measure  the  mean  velocity  in  the  laminar  sublayer 
ai;u  the  s iatic - pre s sur e gradient  »n  the  direction  of  flow. 

Several  pressure  measuring  devices  with  greater 
sensitivities  have  been  reported  in  literature,  e.  g.  reference  3L 
But  generally  they  are  difficult  to  construct.  Laufer  in  re- 
ference 6 stated  that  he  obtained  a sensitivity  of  10”  cm.  of 
water  .with  an  inclined-tube  manometer  using  Benzol,  and  a cathe- 
tometer  to  measure  the  meniscus  displacement. 

The  inclined  manometer  is  attractive  because  of  its 
simplicity.  The  measurement  is  reduced  to  a simple  deter- 
mination of  length  to  normal  engineering  accuracy.  Its 
accuracy  is  limited  by  complications  associated  with  the 
capillary  force  on  the  meniscus. 

A cursory  literature  search  did  not  uncover  an  analysis 
of  either  the  static  or  dynamic  response  of  the  reservoir-type 
inclined  manometer.  Therefore  a preliminary  analysis  was  made. 
The  result  was  that  the  static  problem  appeared  to  hinge  upon  the 
reproducibility  of  the  capillary  rise,  e.  g.  •,  surface  cleanliness, 
temperature,  etc. 

A piiot-model  manometer  was  assembled  and  it  was  found 
that  the  above  difficulties  were  tractable.  Therefore  a final  mano- 
meter was  constructed  and  this  is  shown  in  figure  34. 
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The  basis  of  the  liquid-type  manometers  is  the  pres- 
sure-head relation 


= ft,  - h)  h 
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where  one  column  of  fluid  has  the  density  Jf,-  for  a height 
A , the  other  has  the  density  for  the  same  height,  and 
over  the  remaining  height  of  the  two  columns  from  the 
point  of  inter-conuection  to  the  (assumed)  common  elevation 
where  /If*  i?  measured  the  density  is  the  same  in  each 
column.  Obviously  this  could  be  generalized  to  include 
variable  densities  in  the  columns  if  one  of  the  fluids  is  com- 
pressible , such  as  a gas.  And,  indeed  these  cf.c c iu  may 
not  be  negligible,  but  are  generally  easily  computed  and  will 
not  be  discussed  here. 

The  vertical  (/•  tube  manometer  is  the  simplest  but 
H is  difficult  to  measure  to  the  accuracy  desired.  The 
differential  manometer  makes  the  term  { S,  ~ St)  small  to  en- 
large H but  difficulties  have  been  encountered  with  the 
meniscus  interface.  The  inclined  manometer  attains  its 
accuracy  by  making  the  actual  length  measured  many  times 
larger  than  h . It  appeared  to  be  the  most  probable  for 
success,  so  it  was  selected. 

Distilled  water  was  used  in  this  experiment  because 
of  the  ease  of  maintaining  constant  density.  Alcohol,  for  ex- 
ample, is  significantly  hydroscopic.  Furthermore  the  thermal 
expansion  of  water  is  about  one-fifth  that  of  other  liquids.  The 
disadvantage  of  water  is  its  high  surface  tension.  This  was 
reduced  from  73  to  30  dynes  per  cm.  , a value  typical  of  many 
liquids,  by  adding  a wetting  agent  (The  Mc-riam  Instrument 
Company,  D-2930)  which  does  not  change  the  density. 

The  chief  question  concerning  the  inclined  manometer 
is  whether  the  capillary  rise  is  repeatable  and  constant  over 
the  length  of  the  tube.  Therefore  this  will  be  pursued.  The 
following  variables  determine  the  capillary  rise:  the  diameter 
and  inclination  of  the  tube,  the  liquid,  and  the  angle  of  contact 
between  the  liquid  and  the  tube.  We  would  now  like  to  find  how 
each  of  these  affects  the  capillary  rise. 
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Intuitively  it  see  me  that  there  is  an  optimum  menis- 
cus shape  (the  rise  of  each  point  on  the  surface  beyond  the 
zero-capillarity  position)  to  minimize  capillary  effects,  i.  e. 
to  have  a minimum  sensitivity  to  surface  wetness,  slight 
changes  in  angle  or  surface  tension,  or  whatever  else  may 
change  the  position  of  the  meniscus  minimum.  For  example, 
as  the  tube  inclination  is  increased  the  meniscus  spreads  out 
further  and  consequently  it  would  appear  to  be  more  sensitive 
to  the  way  the  contact  is  made  at  the  glass.  But  at  the  same 
time  the  sensitivity  to  pressure,  which  is  what  we  are  looking 
for,  increases.  The  ideal  procedure  would  be  to  obtain  an 
analytical  expression  for  the  meniscus  shape  then  proceed  to 
optimize.  Unfortunately  this  cannot  be  done  by  a membrane 
analysis  because  the  angle  of  contact  between  most  liquids, 
including  water,  and  clean,  wetted  glass  in  the  presence  of  air 
is  reported  to  be  O degrees,  reference  32  . In  the  present 
case  the  meniscus  shape  and  tube  inclination  were  selected 
by  inclining  several  tubes  at  about  a 20:1  slope  and  then  taking 
the  one  with  the  largest  tube  diameter  that  did  not  have  an  un- 
duly long  meniscus.  Obviously  some  study  could  be  given  to 
selecting  an  optimum  meniscus  shape  and  tube  inclination. 


However,  having  selected  a meniscus  shape  and  tube 
inclination  we  now  wish  to  keep  them  fixed  while  we  examine 
the  effect  of  the  other  variables.  We  can  do  this  by  obtaining 
a similarity  parameter.  The  forces  involved  are  the  gravi- 
tational force  and  capillary  force,  no  we  shall  call  their  ratio 
M.  This  parameter  was  not  found  in  the  literature  but  it  can 
be  obtained  from  the  ratio  of  the  Weber  number  to  the  Froude 
number 


fW f 
LFr 


Weber  //o. 
Froude  Ho. 


l7=  , 

W jh 

J r*  v 


111  = M 


where  t is  the  tube  inside  radius,  ^ is  the  local  gravitational 
acceleration,  and  S is  the  surface  tension.  M would  then  be  a 
constant  for  one  meniscus  shape,  one  angle  of  contact  between 
the  meniscus  and  glass,  and  one  tube  inclination.  Henceforth 
we  will  keep  M fixed. 
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As  mentioned,  an  equation  for  the  meniscus  shape 
is  not  available  from  which  we  can  compute  the  position  of  the 
meniscus  minimum.  But  we  can  determine  the  "average" 
height  of  the  meniscus  due  to  capillarity,  /)(  . Define  he  as 

the  product  of  the  volume  of  fluid  raised  by  the  capillary  effect 
times  the  sine  of  the  angle  of  elevation  of  the  tube,^  .divided 
by  the  cross-sectional  area  of  the  tube,  figure  35  , For  a 

verticai  tube  (ije  (j  * 90  degrees)  f\ e can  be  computed  from  a 
balance  of  capillary  forces  and  gravity  for  a O degree  contact 
angle 

(16) 


A 7 mm.  I.  B.  tube  was  chosen  for  the  water -wetting  agent 
combination  in  the  present  case.  In  the  vertical  position  the 
measured  rise  of  the  meniscus  minimum  was  0.  045  in.  By 
equation  16  this  tube  size  and  fluid  would  give  he  - 0.  068  in. 
This  warrants  the  use  of  equation  (16  ) for  the  order-of-magni- 
tude  considerations  that  will  be  made  below,  (The  approxi- 
mation is  actually  somewhat  better  than  appears  since  ht  is 
the  mean  height  somewhere  between  the  meniscus  minimum  and 
maximum.  ) 


Before  going  through  the  analysis  with  equation  (16  ) 
let  us  show  how  this  equation  can  also  be  used  to  obtain  the  cap- 
illary rise  in  an  inclined  tube.  This  will  be  carried  out  by  refer- 
ence to  figure  36  . The  top  sketch  represents  a sectional  view 

through  the  plane  containing  the  axis  of  the  tube.  An  element  of  the 
glass  surface  containing  the  line  of  contact  of  the  meniscus  is  shown  in 
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the  lower  diagram.  Since  the  angle  of  contact  is  0*  , the 
resultant  force,  dT,  on  the  element  of  the  meniacua  peri- 
meter, da,  ie  in  the  plane  of  the  glass  and  proportional 
to  ds. 


dn  sds 


The  component  in  the  direction  of  the  axis,  dT  , is 

z 


dT* 


- fofo 

ds 


dT 


in  polar-cylindrical  coordinates  about  the  tube  axis.  The 
total  force  in  the  axial  direction,  Tz,  ia 

T.  - f S d?dt 

‘ j ds 
around 

t*e  men: scui 
perimeter 


Therefore,  Tz  ia  independent  of  the  tube  inclination.  The 
capillary  force  balances  the  gravity  force  so  that  the  volume 
raised,  Vc  , is  related  by  the  angle  of  elevation  ^ , 

Tjf  = TVe  ’ sin  (3 
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and  by  definition  of 
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//  / 


1L 

fie 


Sin  (3 


then  combining 


T?  = irr'k  r 


so  hc  is  independent  of  the  tube  slope  (as  was  h ).  Of  course, 
the  rise  of  the  meniscus  minimum  in  a vertical  tube  is  some- 
what less  than  hc the  two  becoming  equal  as  r approaches 

O.  But  within  this  difference,  the  rise  of  the  meniscus  min- 
imum in  an  inclined  tube  can  be  approximated  by  the  rise  in  a 
vertical  tube.  This  will  be  done  in  the  following  analysis. 

Let  us  now  select  the  liquid.  Introducing  M into 
equation  ( 16  ; 

r 2 Lmr]1 


M is  fixed,  so  one  chooses  a fluid  with  the  smallest  value:  of 
to  give  the  least  capillary  rise.  A survey  of  a large 

number  of  liquids  revealed  that  Methanol  has  the  lowest  value, 
and  this  is  only  four  percent  below  the  water  - wetting  agent 
combination. 

The  uniformity  requirement  for  the  tube  diameter 
can  be  established  by  putting  equation  (16  ) in  the  form 


_ _ 2$  a _ Z. 
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One  sees  that  all  liquids  under  similar  (which  could  mean  op- 
timum) conditions  have  the  same  sensitivity  to  absolute  varia- 
tions in  radius.  With  Lhe  present  combination  a A r = 0.0003  in 

gives  Ahc  = ft.  0001  in.  Normal  precision  tube  tolerance  is 
->■0.0004  in.  from  the  nominal  size,  and  because  of  the  method  of 
manufacture  the  uniformity  along  one  tube  is  significantly  better. 
Similarly,  the  reservoir  accuracy  was  established  and  turned  out 
to  be  easily  fulfilled. 


Temperature,  0.  enters  in  three  principal  ways: 

(1)  volumetric  change  of  the  container,  (2)  expansion  of  the  liquid, 
and  (3)  variation  in  capillarity.  These  will  now  he  discussed. 

(1)  Consider  a thin-wall  cylindrical  glass  reservoir  of 

circular  cross-section  resting  on  its  base.  9 The  volume  of  liquid 
enclosed  equals  the  depth  x (circumference  rAir  . If  the 

liquid  volume  is  constant  then  for  a thin-wall  reservoir  a one  per 
cent  change  in  circumference  corresponds  to  a two  per  cent  change 
in  depth.  The  thermal  coefficient  of  linear  expansion  of  glass  is 
nearly  5 x 10~^/°F  (reference  33  ).  Therefore,  if  the  depth  is 
about  one  inch,  a temperature  rise  of  1®F  corresponds  to  a de- 
crease in  depth  (zero-level)  of  10"^in.  of  water,  which  is  negli- 
gible. The  changes  in  the  glass  tube  affect  the  level  even  less  if 
its  portion  of  the  enclosed  volume  of  liquid  is  significantly  smaller 
than  that  in  the  reservoir,  as  it  was  in  the  present  manometer. 

(2)  The  thermal  coefficient  of  volumetric  expansion  of 
water  is  0.  0001/°F  (reference  33),  (For  most  other  liquids  it  is 
about  five  times  larger.  ) Thus  the  zero  level  in  a one-inch  deep 
reservoir  would  rise  10"^  in.  per  degree  F.  This  iB  the  accuracy 
desired  so  the  temperature  must  be  controlled  within  this  limit. 

The  above  considerations  are  all  connected  with  the  change 
in  zero  level.  These  effects  can  be  reduced  by  periodically  taking 
zero  readings.  This  procedure  has  been  adopted  during  tests.  How- 
ever, there  is  another  effect  due  to  the  thermal  expansion  of  the 
liquid;  namely,  h is  inversely  proportional  to  T in  the  pressure-head 
relation.  This  causes  a percentage  error,  namely  = 0.  0001 /F". 
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(3)  The  temperature  effect  on  the  capillary  rise  is 
computable  by  putting  equation  {2)  in  the  form 

1 Ms  - ±<>J  ..  i M 
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The  first  term  represents  the  change  due  to  surface  tension 

and  the  other  is  the  density  change  within  fixe  capillary-rise 

portion.  For  the  water -wetting  agent  combination  XAJt 

hfQQ 

- 0.  OOll/F0.  A slight  drop  in  the  zero  level  was  observed 
after  turning  on  the  light  illuminating  the  meniscus.  However, 
within  one  hour  this  effect  tended  to  reach  equilibrium.  Actually, 
it  was  very  easy  to  shield  the  manometer  against  all  thermal 
effects  by  enclosing  it  in  a plastic  bag. 

A side  poin'.  could  be  mentioned.  It  is  conceivable  that 
the  manometer  could  he  designed  so  that  the  thermal  expansion 
effect  of  the  liquid  could  be  nullified  by  the  surface  tension  re- 
cession and  expansion  of  the  reservoir. 

The  manometer,  figure  34  was  constructed  as  follows- 
The  inclined  tube  is  7 mm.  I.  D the  reservoir  is  100  mm.  I.  D.  , 
and  the  tube  inclination  is  20:1.  The  microscope  is  mounted  in 
the  vertical  plane  so  that  light  rays  pass  through  normal  to  the 
glass.  The  meniscus  appears  symmetrical  and  nearly  flat  over 
a wide  region,  and  is  sharp  and  clearly  identifiable  to  1/10  of  the 
smallest  reading  desired.  A dial  indicator  graduated  in  10“-in. 
indicates  the  microscope  travel:  0.  002  ir_.  represents  a pressure 

of  10”4  in.  of  water.  An  internal  scale  within  the  microscope 
has  graduations  equivalent  to  a pressure  of  10  "^  in.  of  water  which 
facilitates  averaging  pre  s sure  fluctuations. 

It  was  found  important  in  the  experimental  technique  to  wet 
the  surface  ahead  of  the  meniscus.  (Reference  32  states  that  the 
contact  angle  of  a non-wetted  surface  is  different).  The  tubes  were 
first  thoroughly  cleaned  with  a standard  acid  cleaner.  Then  they 
were  flushed  with  tap  water  and  finally  distilled  water.  Nevertheless, 
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the  readings  could  be  in  error  by  as  much  as  C.  0004  in.  of 
water  unless  the  meniscus  was  displaced  far  enough  ahead 
to  wet  the  gla.ss.  With  the  wetting  technique,  the  meniscus 
would  return  to  the  same  position  with  the  probable  error 
of  a small  fraction  of  10"^  in.  of  water.  In  this  connection, 
one  must  be  careful  to  avoid  foreign  matter  since  it  generally 
tends  to  accumulate  at  the  surface. 


Evaporation  effects  can  enter.  The  zero  level  with 

A 

the  reservoir  open  was  observed  to  drop  by  10"*  in.  in  a 
five-minute  period.  Therefore,  a constriction  was  found 
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vent  evaporation  from  changing  the  zero  level,  and  large 
enough  so  the  vapor  pressure  does  not  differ  appreciably  from 
the  ambient.  For  example,  the  zero  level  consistently  rose 
0.0002  in.  when  the  reservoir  connection  was  switched  from  a 
1/4  in.  I.  D.  tube  open  to  the  atmosphere  to  an  identical  tube 
with  a 0.  020  in.  diameter  hole  lorming  a constriction  to  the 
atmosphere.  This  is  the  pressure  difference  due  to  the  flow 
of  vapor.  A slightly  larger  opening  gives  no  effect. 


Three  calibration  tests  were  made: 


(1)  Zero  reading  - A number  of  zero  readings  were 
taken  and  the  spread  was  always  within  10”^  in.  of  water  so 
long  as  the  fluid  was  sufficiently  clean.  Following  several 
weeks  use,  the  fluid  would  become  dirty  which  always  mani- 
fested itself  by  not  duplicating  zero  readings. 

(2)  Hysteresis  - Small  quantities  of  liquid  were  added 
in  the  reservoir  and  the  rise  was  measured  in  the  tube.  Then 
the  quantities  were  successively  removed.  No  hystersis  effect 
was  observable. 

(3)  Accuracy  - The  final  calibration  was  made  by  adding 
measured  volumes  of  liquid  in  the  reservoir  by  a precision 
burette.  The  probable  error  over  a range  equivalent  to  0.  2 in. 
of  water  pressure  was  0=  00007  in.  of  water. 
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APPENDIX  III 

CORRECTION  OF  THE  HOT-WIRE  SENSITIVITY  DUE 
TO  THE  TEMPERATURE  DEPENDENCE  OF  KING'S 

CONSTANT,  A 


As  mentioned  on  page  23,  the  so-called  constant.  A,  in 
King's  equation  (ref.  22)  is  3. c t u 3. 1 lv  a function  of  the  wire  temper- 
ature.  This  was  originally  observed  by  King,  but  evidently 
has  been  overlooked  in  a nember  of  intervening  investigations. 
Therefore  a correction  must  be  made  when  the  hot-wire  is 
calibrated  at  constant  resistance  (constant  temperature)  and  used 
at  constant  current.  In  the  present  experiment  this  correction 
(a  monotonic  function  of  overheat  ratio)  increased  the  computed 
values  of  ail  U.  by  16  to  24  percent.  Therefore  the  correction 

7 T 

will  be  derived  below. 

King's  equation  for  the  heat  loss,  I^R;  of  an  infinite  circular 
cylinder  with  axis  normal  to  the  stream  direction  is 


ISH  = [a(-O)  * t%)  07) 


where  U is  the  air  speed,  constant  and  low,  "$■  is  the  wire  temper- 
ature , and  is  the  ambient  air  temperature.  A(t9")  is  the 
term  under  examination,  but  for  some  operating  conditions  the 
temperature  dependence  of  R is  equally  important  so  it  will  be  in- 
cluded. King's  result  that  B is  only  slightly  temperature  sensitive 
has  been  verified,  so  it  will  be  assumed  constant. 

To  apply  equation  (17)  to  turbulence  measurement  the  usual 
assumption  of  a quasi-equilibrium  state  of  the  flow  field  will  be  admitted. 
The  justification  will  not  be  undertaken  here.  So  for  constant- 
currenL  operation  of  the  hot-wire  in  a constant-^  stream  the 
functional  form  of  equation  (17)  is 

R - <R[A(-&)  , t 9 , U] 
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-&  is  independently  related  to  R by  the  thermal  resistivity,  so 
neglecting  second-order  and  higher  terms  in  a Taylor's  expansion 
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The  term  £ embodies  the  correction,  and  usually  has  been 
neglected  without  explicit  mention. 

i(Q  ^ ft  J ^ 

The  terms  ?■  * . t-  • , and  ",  are  easily  obtained 

lu  1a  J# 

from  equation  (17)  by  differentiation.  For  A(tJ^)  we  will  use  the 
relation  defined  by  King 


A = Ae  ft  + c (■&  - -<.%}] 


which  gives 
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The  evaluation  of  c will  be  discussed  later.  The  remaining 
tp'm.  dR  . is  obtained  from  the  customary  expansion 


R = Rf[i  * + 


v/here  is  the  reference  temperature  lor  R,  usually  different 

w.e  • 

Substituting  these  expressions  into  equation  (18)  gives 
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(20) 


We  will  now  evaluate  A/D  from  a calibration  of  the  hot-wire  at  a 
constant  resistance,  figure  37.  The  value  of  the  resistance  in 
curve  I corresponds  to  the  mean  operating  resistance.  This  curve 
will  be  used  for  the  present  calculation.  We  will  adopt  several 
convenient  definitions  from  reference  21. 
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where  vu,  corresponds  to  the  abc.issa-intercept  of  curve  I. 
Equation  (20)  becomes 
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Using  the  following  definitions  of  "overheats" 


R-  Rf 
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and  Ohm's  law,  we  obtain  for  constant  current 


AC  __  _ qS  2 AU 
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(21) 
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whore  E is  the  voltage  drop  across  the  hot-wire  and 
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tiy  the  following  approximations, 
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The  expression  for  reduces  to 
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Equations  (Zl)  ^with  (22) J expresses  the  hot-wire  response  at 
a constant  current  to  velocity  fluctuations.  Equatio  i (21)  is 
identical  to  equation  (12)  of  reference  ?l  when  cd  - 0. 

The  effect  of  upon  the  hot-wire  response  is  shown  in 
figure  37.  (1)  For  A - constant  the  straight  lines  I,  II,  and  III 

would  represent  three  typical  cons tant- resistance  responses 
*11  > *1  > Rjjl  of  a hot-wire  following  King's  equation.  The 

constant -cur  rent  response  of  the  same  hot  wire  would  be  repre- 
sented by  the  dashed  line  , h,  between  curves  II  and  III. 

(2)  In  the  actual  case  where  A - A{q9*  ) the  straight  lines  I,  IT', 
and  III1  represent  the  con  stant  - re  s i stance  response,  and  the  constant- 
current  operation  would  follow  the  dashed  line,  j , between  II'  and 
III1.  Obviously  since  R T j , and  Rjjj  R ITT • t^e  hot-wire 

signal  is  reduced.  Or,  in  other  words,  for  a given  signal  from  a 
hot  wire  in  a fluctuating  velocity  field  the  true  velocity  fluctuation 
is  greater  than  would  be  computed  by  assuming  J = 0. 

King  found  empirically  that  with  a wire  perpendicular  to  the 
flow,  c = U. Guild  Cent,  degrees  for  air  based  upon  = 17 
deg.  Cent.  This  value-  was  verified  in  rhe  present  experiment  with 
a 0.00015  in.  dia.  platinum  wire  0.040  in.  long  in  a non-turbulent 
air  stream  at  15  to  35  feet  per  second.  The  wire  was  oriented 
in  two  positions  relative  to  the  flow,  (i)  normal,  and  (2)  inclined  50 
degrees  from  the  normal.  The  procedure  was  to  calibrate  the  wire 
at.  a‘w  = 0.4.  0.7,  and  i.O  in  each  of  the  two  positions.  The  value  of 
A %vas  computed  by  assuming  B was  independent  of  . The 
comparison  of  c to  King's  value  is  made  in  figure  38.  Case  l is 
in  agreement  with  King's  value.  The  trend  in  case  2 may  represent 
experimental  error.  This  warrants  further  investigation. 

The  values  of  and  (3  for  platinum  and  tungsten  were  taken 
from  reference  32  (Vol,  VI,  page  136)  and  are  tabulated  beiuw  ir. 
the  form  in  which  they  appear  in  equation  (22). 
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The  following  remarks  pertain  to  equation  (22). 


1. 


oj. 

c 


dominates  the  denominator  of  the  first  tern!. 


This  is  the  A(t9-)  effect.  Therefore  the  values  of  are  general- 
ly larger  for  platinum  than  for  tungsten.  { 2 is  a function  of 
aerodynamic  properties  and  wire  geometry,  and  is  independent 
of  the  wire  material.  ) 


- a.  . . .... 

is  negative  for  platinum  and  tlieretore  increases  j . 
whereas  for  tungsten  this  term  is  positive  and  tends  to  reduce  <j>  . 

In  the  Dresent  exDeriment  contributed  about  20  Dcrcent  to 

Of*  t 

the  value  of  cb  , or  4 percent  to  the  value  of  U. 

, ^ U 

3.  For  values  of  Q ^ < O.S"  , % is  roughly  proportional 
to  the  overheat,  . Of  course  the  signal  is  also  roughly 

proportional  to  c so  it  is  not  necessarily  feasible  to  reduce  Qv 
in  order  to  reduce  jp  . 


In  figure  39  the  <j£  correction  has  been  plotted  for  both 
platinum  and  tungsten  as  a function  of  the  velocity  term,  2 , with 
Clur  as  a parameter.  The  difference  between  platinum  and 
tungsten  is  evident.  In  the  present  experiment  tungsten  wires 
would  have  reduced  the  values  of  cT  by  a factor  of  2. 


A few  remarks  should  be  made  concerning  the  hot-wire  end 
effect  on  . Equation  (17)  presumes  an  infinitely  long  circular 
cylinder.  The  finite -length  hot  wire  does  not  have  a uniform  temper- 
ature  distribution.  But  simple  considerations  would  indicate 
that  the  average  (spatial)  iff-  corresponding  to  the  average  (spatial) 

R is  probably  satisfactory  for  equation  (17)  with  the  usual  hot-wires. 
But  the  heat  conducted  out  the  end  of  the  wire  to  the  holder,  a length 
effect,  is  net  of  the  form  of  equation  (17).  Presuming  a fixed  ternper- 
ature  of  the  holder  at , and  estimating  the  heat  conducted  out  of 
the  end  of  the  wire,  Q]g  , to  be 


- k (-&  - 

and  taking  first  order  temperature  changes  in  k 

k = ks[i  cK(&  - -&J] 
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it  is  evident  that  the  problem  is  to  compare  c to  c^  • 
Reference  32  (Vol.  V,  p,  221)  gives  the  following  values 
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Since  for  platinum  c^  anu  c are  of  the  same  order,  and  in  addition 
the  end  losses  are  probably  small  for  most  hot-wires,  the  end  losses 
could  be  included  in  A.  For  very  short  hot-wires  this  should  be 
reexamined.  But  for  tungsten  the  correction  is  of  the  opposite  sign. 
Therefore  for  tungsten  the  end  losses  should  be  either  isolated  in 
equation  (17),  or  be  shown  to  be  negligible, 


I 
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APPENDIX  IV 

Remarks  on  the  Calibration  of  the  X-meter 

The  X-meter  was  used  to  measure  v*,  ,/,  uv  and  inciden- 

tally, u'.  It  turns  out  thfit  these  measurements  are  subject  to 
appreciable  instrument  errors  under  ordinary  precisions  of  manu- 
facture and  operation.  In  particular  the  measured  uv  is  very 
sensitive  to  a small  error  in  the  probe  alignment.  Therefore  (1)  an 
analysis  will  be  made  of  the  effect  of  a first-ordcir  difference  be- 
tween the  wires  in  the  sensitivity  to  u and  v (fluctuations)  in  the 
light  of  existing  data.  (2)  A procedure  to  nullify  the  effect  of  the 
differences  will  be  out-lined.  The  '*  black  box“  technique  of  cali- 
brating the  X-meter  described  below  is  similar  to  the  procedure 
described  in  reference  23.  The  present  report  extends  this  work 
by  items  (1)  and(Z)  above. 

The  X-meter,  figure  40,  consists  of  two  fairly-well 
matched  wires  oriented  in  an  “x“  somewhat  symmetrically  about 
the  mean  stream  direction,  i.e.  *Pt.  For  small  fluctuations, 

the  conventional  linearized  hot-wire  analysis  permits  writing  the 
voltage  sensitivity  of  wires  no.  1 and  no.  2 to  fluctuations  u and 
v in  the  form 


~ d f A4.  -f*  h,  V 
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where  the  sensitivity  coefficients  are  in  general  not  equal.  De- 
fine the  differences 


A flu  “ a x~  Q» 
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where  the  non-zero  values  of  A a and  Ah  represent  the  inevitable 
physical  and  geometrical  differences  between  the  two  wires. 
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In  order  to  compute  the  turbulence  terms  from  the  voltage 
measurements  we  will  assume  4«  a and  b.nnd  neglect 

second-order  terms  (c.g.  ),  obtaining 
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In  order  to  compute  the  turbulence  quantities  we  must  determine 
a,  b,  A a,  and  .fib.  We  will  discuss  the  first  two  later.  Two  alter- 
native procedures  are  available  in  connection  with  a a and  A b.. 
(l)Aa  and  ab  could  be  measured  at  whet  ever  value  they  happen 
to  have  and  equations  (24),  (2S),  and(26)  would  then  be  solved  simul- 
taneously for  the  unknowns  u,  v and  uv,  (2)  The  probe  could 
be  reoriented,  etc.  to  make  a a = 0 and  ab  = 0 and  then  the 

equations  would  each  have  only  one  unknown.  Our  experience  has 
been  that  it  is  more  expedient  to  use  method  (2),  without  sacrificing 
accuracy,  than  to  measure  the  original  values.  Therefore  the 
second  procedure  has  been  used,  and  we  will  focus  our  attention  on 
examining  the  effect  of  the  error  involved  in  ascertaining  that  A a 
and  A b are  indeed  zero. 

The  three  conditions  that  the  final  X- meter  will  satisfy  are: 


A CL  = O 


U)  ,&  b = o 


(3)  M,  ~ Mc 


where  M , is  the  time  constant  of  wire  no.  1. 
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Tn  order  to  satisfy  condition  (3)  wt  will  overheat  the  wires  equal- 
ly, i.e.  = ct ^ . The  overheat  condition  does  not  have  to  be 

satisfied  »-  exactly  as  conditions  (1)  and(2).  We  have  three  adjust- 
able quantities  available  to  fulfill  these  requirements: 

(1)  the  probe  orientation,  «", 

(2)  a voltage  divider  across  the  output  from  one  v/ire 

(3)  the  current  into  each  wire  (by  a resistance  in  series 
with  one  not  wire  if  it  is  in  a bridge,) 

Arbitrarily  wc  will  place  the  voltage  divider  across  hot  wire  no,  2 
so  that  its  measured  voltage  is  tetj  . Henceforth  the  notation  will 
be  changed  so  that  A a = ka  - a and  4b  = kb  - b . 


We  will  now  estimate  the  effect  of  an  angle  change  on  the 
a*s  and  bs.  The  assumption  will  be  made  that  if  a stream  approaches 
a wire  obliquely  it  is  the  component  normal  to  the  wire  that  does 
the  cooling.  This  is  not  accurately  borne  out  by  experiment,  but  it 
is  close  enough  lor  a good  estimate.  We  will  not  use  it  to  evaluate 
the  coefficients  , but  only  to  estimate  the  order  of  their  change. 

King's  equation  for  this  inclined  wire,  figure  40,  becomes 
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where  y » \fy  , 5 (/cos  , y„.  - , and  the  second- 

order  resistance  coefficient  is  neglected  so  and  in- 
clude only  the  first-order  resistance  coefficient  oc  (in  contrast 
to  Appendix  Til  ).  (Note:  the  subsciipt  i will  be  suppressed 
through  the  remainder  of  this  paragraph)  After  several  algebraic 
steps  equation  (28)  becomes 
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where  E is  the  voltage  across  the  wire.  Define 
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Tin'll  the  coefficients  a and  b (for  the  single  wire)  are  for 
constant  current  operation 
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The  probe  is  rotated  by  a * to  a new  position.  At  the  new  position 
is  made  the  same  by  changing  the  current  (although,  of  course, 
once  in  the  new  position  the  X-meter  operates  at  constant  current). 
Express  E as 
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The  indicated  steps  give 
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Note  that  0<  i < 1 and  -p-=  = ± 1 so 
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At  the  wall  where  U—+  v)  2 — ► 0;  as  U-*  «e  , f 1. 

Roughly  b is  3 to  7 times  as  sensitive  to  the  probe  orientation 
as  a.  At  45*  equations  (29)  give  a = b. 

We  now  have  sufficient  equations  to  evaluate  the  effect  of 
an  angle  change  of  the  probe  on  the  turbulence  measurement,  For 
simplicity  the  two  wires  will  be  oriented  at  , - t*  = 45  degrees 

and  <f  = 0 degrees.  (However,  the  equations  are  general  and 

t tie  extension  is  obvious.)  The  angle  change,  * f , affects  the  tur- 
bulence terms  through  A a and  A b in  equations  (24),  (25),  and 
(26).  From  equations  (30#) 
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and  within  this  approximation 
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= 1 and  "3^  = ^8ure  4-0.  Therefore  near  the 
wall  a 1 degree  error  in  Y corresponds  to  a 1.7  percent  error 
in  and  a 5 percent  error  in  . In  this  region  the  effect  on 

u£  is  small  since  -^tX-  *■<  1;  on  v^*"  also  small  as  pi  1;  hut 

— r *•  L .fj 

•*<  - -»  ~>  n ir'..  „ i=. „ _ - , - v,  i . In  in  = n{j  . 


> 20  near  the  laminar  sublayer  (figure  30  ana  reference  6) 
one  degree  t^r^or  would  give  an  error  in  uv  of ' 15  percent. 

- i .-i  a >0.*  / 1 _L  i r i.i  ; r 


so  a 

At  y * = 100,  -4=^-  =l  6 and 

a one  degree  error  would  be  approximately  2 percent. 


1.5,  the  error  in  uv  for 


A method  of  adjusting  the  probe  io  satisfy  equations  (27) 
and  measure  a and  b will  nov;  be  presented.  Two  general  method 
(of  unequal  accuracy)  are  obviously  available 

(1)  ~ Cosine-iaw“  method,  for  approximate  measurements: 
each  wire  is  calibrated  separately  normal  to  the  air  stream  in 
the  usual  constant-resistance  way.  Then,  by  finding  the  = 0 
positions  and  using  equation  (28)  the  terms  in  equations  (29)  are 
computed,  giving  aVtand  and  thus  determining  the  position 
which  satisfies  equations  (27).  In  addition  to  inaccuracy  of  the 
“Cosine  law“f  equation(28),  a weakness  of  this  method  *rises  in 
the  error  of  determining  the  ft,x  - G position  for  each  wire.  (The 
effect  of  wire  differences  has  been  analyzed  by  reference  3 where 
the  wire  sensitivities  are  determined  by  the  “Cosine-law“  method.) 

(2)  “ Black-box“  technique  for  more  accurate  measure- 
ments: Pragmatically  consider  the  X-meter  as  an  electrical 
device  which  simply  responds  to  U and  V in  some  unknown  way. 
(u,  v,  and  e are  increments).  Keeping  only  the  first-order  terms 
5n  a Taylor  series  expansion,  we  obtain  equation  (23),  as  before. 
But  this  time  we  will  find  the  sensitivity  coefficients  empirically, 
by  installing  the  X-meter  in  a stream  where  u and  v are  varied 
independently.  During  the  calibration  we  measure  only  C,  + 

and  . 


To  first  satisfy  equations  (27)  the  following  procedure  is 
employed.  Place  the  probe  in  a calibration  duct  at  roughly  the 
“neutral  angle"  (where  & a = a b 0)  with  the  correct  a^. 
Note  that 
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The  obvious  procedure  is  to  keep  the  probe  fixed  and  vary  the  air 
speed  (at  “constant  current11  ),  then  to  fix  the  air  speed  and  rotate 
the  probe.  (Actually  it  will  be  apparent  that  the  process  is  hastened 
by  measuring  A b first.)  Typical  curves  are  snown  in  figure  41. 

To  satisfy  equations  (2?)  we  can  change  if  , and  utilize  the 
voltage  divider  measuring  ke4  (or  ke(  ) where  the  first  measure- 
ment was  then  at  k = 1.  The  value  c,f  a^,_  is  the  same.  Indicate 
the  second  measurements  by  a prime  {/)  Using  - Taylor  scries 
expansion  to  fii*>t  order  terms  only 


ao.  5 ~ -b  (V~  r)  + d x(h'~  it.) 

A b 5 it  jat  - - k i-  Jt  (~^  ) , (*  ~Jr)  f ^ 

tur 
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This  is  solved  lor  AX'  and  with  Aa'  = A b = 0. 
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£ Simulating  operation,  i.e.  without  adjusting  electronic 
equipment. 
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The  coefficients  of  the  sec  ond  terms  on  the  right  hand  side  are 
evaluated  by  equations  (30).  Generaily  the  approximate  values 
of  equations  (30*)  will  suffice.  Equations  (30  ) give  for  the  X- 
meter  of  figure  40 

o ~ (kex-ot)  * j AY  + ~— 


c s (k.  a,  ) - (ft  + )( - J A V h. 

The  new  values,  (k  a 2 - a,)  and  (k  - b4 ) , should  be  mea- 
sured by  repeating  the  foregoing  procedure,  to  verify  that  equa- 
tions (27}  are  indeed  satisfied,  and  the  process  could  be  repeated 
again  although  generally  it  was  not  found  necessary. 

The  nexc  step  in  iue'diou  (2)  is  io  measure  a and  b.  ihe 
method  is  obvious  following  the  procedure  for  a a and  a b since 


= T TZ  (k  * C'  ) 
h a"i  Jv 


The  accuracy  in  determining  a and  b is  very  good,  and  there  car> 
be  little  doubt  as  to  the  correctness  since  the  probe  is  calibrated  in 
the  manner  in  which  it  is  used.  In  principle  this  calibration  should 
be  repeated  over  the  range  of  mean  velocities  required.  However 
U enters  through  i!  and  if*  of  equations  (29).  Calibrating  the  hot- 
wire at  constant  R and  treating  the  two  wires  as  one  gives  y0  , 
and  thus  ? and  ^ can  be  computed. 

In  practice  it  was  noted  that  the  calibration  is  limited  in 
accuracy  by  the  determination  of  A.  b,  The  velocity  is  fixed,  the 
probe  rotated,  and  e,  + k ee  is  measured.  But  since  the  coeffi- 
cient a is  relatively  large,  e , + ke,  is  sensitive  to  small  varia- 
tions in  the  air  speed  of  the  calibration  stream,  'f’his  is  not  a 
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weakness  singular  to  method  (2)  but  also  occurs  in  method  (1)  in  the 
determination  of  = 0 and  A and  B in  King*  s equation.  This 

kind  of  difficulty  is  less  troublesome  in  the  measurement  of  Ct  a 
since  the  direction  oi  tiie  calibr aiiuu  sli  earn  is  generally  more  con- 
stant. The  difficulty  in  determining  Ab  is  indicated  in  figure  41. 

The  spread  of  the  points  indicates  perhaps  an  uncertainty  in  deter- 
mining the  “neutral  position**  of  about  ±.  1-2  degrees. 

Once  a value  of  i'  had  been  selected  the  X --meter  was 
oriented  in  the  wind  tunnel  by  measuring  ke.  - e, , see  figure  41, 
and  rotating  the  probe  to  the  appropriate  value.  This  could  be  achieved 
within  ± 1/2  degree.  Therefore  it  is  evident  that  the  problem  in  the 
present  case  was  not  resetting  the  probe  to  the  same  angle,  but  find- 
ing out  what  angle  is  the  correct  one  for  satisfying  equation  (27). 

Incidentiy,  in  a boundary-layer  type  flow,  if  the  probe  is 
simply  translated  "laterally**  a realignment  may  be  necessary  near 
the  surface  since  the  streamlines  are  not  parallel. 
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APPENDIX  V 


EQUATIONS  OF  MOTION  IN  CYLINDRICAL- POLAR  COORDINATE 
FOR  FLOW  IN  A T WO  - dtk  rF.  k>  sT  ON  A T . CONVERGING  OR  DIVERGIN’ 

CHANNEL 

Intr  ori  ucti  on 

For  the  steady-state  two-dimensional  laminar  flow  be- 
tween conve  rging  or  diverging  plane  walls  a general  solution  has 
been  obtained  by  Rosenhead  * by  applying  the  restriction  of  radial 
flow.  Rosenhead  discussed  in  detail  the  effect  of  the  Reynolds  num- 
ber and  the  angle  between  the  walls  lor  both  the  inflow  and  outflow 
cases . 


In  view  of  the  success  achieved  for  the  laminar  case  an 
analysis  was  undertaken  at  the  outset  of  the  present  turbulent  flow 
investigation  to  learn  whether  the  restriction  cf  radial  flow  applied 
to  the  mean  motion  might  yield  additional  information,  perhaps  per- 
mitting one  integration  of  the  resulting  ordinary  differential  equation. 
Since  purely  radial  mean  flow  is  presumably  possible  only  with  full 
dynamical  similarity,  it  is  implied  that  all  appropriately  dimension- 
less statistical  properties  of  the  turbulence  must  also  be  functions 
of  the  angle  only.  Unfortunately  the  final  equation  contained  four  un- 
knowns in  the  general  case  which  could  still  be  reduced  only  to  two 
unknowns  by  the  additional  assumption  of  low  turbulence  level.  In 
the  absence  of  an  independent  relation  this  approach  was  temporarily 
abandoned . 

The  next  attempt  was  through  the  integral  relations.  The 
differential  equations  in  cylindrical-polar  coordinates  were  integrated 
from  the  centerline  to  the  wall  along  a constant  radius  arc  similar 
to  the  Cartesian  case  (p.  8).  The  final  equation  was  not  different 
from  the  Cartesian  integral  equations  for  a small  angle  so  this  was 
not  pursued  further. 

It  was  only  near  the  end  of  the  present  investigation  that  a 
pragmatic  advantage  of  cylindrical- polar  coordinates  in  the  diver- 
gent channel  was  realized:  The  Cartesian  coordinates  had  been 
oriented  on  the  centerline  of  the  channel  in  order  to  simplify  the  in- 
tegral equations  (by  incorporating  symmetry).  But  the  disadvantage 
was  that  at  the  wall  these  coordinates  are  not  parallel  to  the  surface,  so 

* Rosenhead,  L.:  The  Steady  Two-dimensional  Radial  Flow  of  Viscous 

Fluid  Between  Two  Inclined  Plane  Walls.  Proc.  Royal  Soc.,  Ser,  A., 
v.  175 . 1940,  p.  436-67. 
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of  turbulent  energy  terms  (footnote,  p.  7,  and  p,  29).  Cylindrical- 
polar  coordinates  eliminate  this  problem  since  they  are  aligned 
with  both  the  centerline  as  well  as  the  wall.  There  was  insufficient 
time  available  when  this  was  discovered  to  rework  the  whole  paper. 
Therefore  the  differential  equation  cylindrical-polar  coordinates 
are  presented  in  this  appendix.  In  addition  the  ordinary  differential 
equation  for  radiai  flow  described  above  will  be  deduced. 


liquations  of  Motion  in  Cyiindrical-Poler  Coordinates 

Consider  the  cylindrical-polar  coordimites  r,  u»  , and  z 
where  w = 0 is  the  centerline  between  the  walls  and  r ~ 0 is 
at  their  intersection;  U,  V,  and  W are  the  instantaneous  velocities 
in  the  respective  directions.  Again  the  time  (or  ensemble)  average 
will  be  i idicated  by  a bar,  etc.,  so 


u f i.j&)  = u (*,  « ) *■  ( * , **, 

etc. 


PfJt'Uu,  = P (A,  o>  i ) + (Jt,  w,  it  , A ) 


The  Reynolds  equations  and  turbulent  kinetic -energy  equations  in 
cylindrical-polar  coordinates  have  been  written  cat  independently 
by  Wislicenus  and  Yeh  (ref.  34)  so  they  will  be  presented  here  merely 
in  their  final  form. 
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The  continuity  equation  (31),  Reynolds  equations  (32),  (33), 
and  (34)  end  turbulent  kinetic-energy  equation  (35)  are 
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The  first  three  groups  of  terms  in  equation  (35)  were  identified 
essentially  as  in  reference  34:  Group  I represents  the  production 

of  turbulent  energy  from  the  mean  motion,  where  the  last  two  terms 
are  associated  with  the  turbulent  centrifugal  and  corioilr 
Group  IT  represents  the  convectio..  of  turbulent  energy  by  the  mean 
flow;  and  Group  III  represents  the  convection  of  energy  by  the  tur- 
bulent motion  plus  the  transfer  of  energy  by  the  work  of  fluctuating 
pressure  gradients.  The  following  identification  is  given  to  the  re- 
maining terms:  Group  IV  represents  the  diffusion  of  turbulent  en- 

ergy by  viscous  forces  (i.e.  viscous  v/ork  terms);  and  Group  V 
represents  the  dissipation  of  turbulent  energy  to  heat.  Group  V 
was  identified  independently  following  the  method  of  Goldstein,  et  al.  * 
extended  to  the  turbuient  motion. 


Radial  Flow 


The  mean  flow  will  now  be  assumed  to  be  radial,  and  all 
mean  values  of  velocity  terms  will  be  assumed  to  be  two-dimensional. 
Specifically 


V = \A/  — — 


. c 


where  f is  any  velocity  term 


o x. 


The  continuity  equation  (31)  reduces  to 
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It  is  evident  that 
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* Goldstein,  S.,  editor:  Modern  Developments  in  Fluid  Dynamics, 
Vol.  T,  First  Edit.,  1938, ‘p.  98. 
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These  are  equivalent  to  the  assumption  ox  longitudinal  similarity 
for  turbulent  terms  in  the  C 
-5^7  , etc.  are  functions  o( 

(34)  reduce  respectively  to 


for  turbulent  terms  in  the  Cartesian  case  (p.  11),  since  4^*  , 

■^7  , etc.  are  functions  of  uj  only.  Equations  (32),  (33),  and 
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Differentiating  equation  (40)  with  respect  to  z and  equating  with 
equation  (39)  gives 
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Differer.tiating  equation  (40)  with  respect  to  A we  obtain 
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Equating  with  equation  (37)  gives 
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and  since  the  rieht-hand  side  is  a function  of  Uf  alone  define 
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fa  a constant 


The  constant,  k.  can  be  evaluated  from  th 
at  the  wall  by  equation  (41)  since 


e radial  pressure  gradient 
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where  i*  is  again  the  half-angle  between  the  walls. 
Therefore 
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Equation  (42)  is 
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Equation  (43)  is  the  primary  differential  equation  for  the  flow  field. 
Its  counterpart  in  the  Cartesian  coordinates  is  equation  (6),  except 
that  equation  (43)  applies  without  restriction  on  ot  or  turbulence 
level.  Unfortunately  there  are  four  unknowns.  The  corresponding 
equation  for  parallel -wall  channel  flow  (equation  (3a)  of  reference  5) 
can  be  integrated  once  directly. 


Equation  (43)  reduces  to  two  unknowns  where  the  turbulence 
level  is  sufficiently  small  that 


-T2.1 


in  which  case 


(44) 


= 0.  Since 

two  unknowns  remain  in  the  equation  an  independent  relation  is  re- 
quired for  its  solution.  This  was  not  available  so  the  analysis  was 
terminated . 

The  turbulent  energy  equation  (35)  is  not  significantly- 
simplified  in  radial  flow.  The  eleven  production  terms  reduce  to 
three,  the  three  convection  terms  to  one,  but  only  two  of  the  remaining 
terms  in  III,  IV,  and  V are  eliminated.  Therefore  the  turbulent 
energy  equation  will  not  be  presented. 


IX, 


ft 


-TL 


rr  ( -a,*  - .a'  ) = h 


Equation  (44)  is  identical  to  the  laminar  case  when 
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Diagram  of  two-dimensional,  channel* 


Figary  4.-  Block  diagram  of  turbulence  measuring  equipment. 
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Figure  6.-  Mean-velocity  distribution  in  1 in.  channel. 
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Figare  5.-  Kean- velocity  distribution  near  wall 
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Figure  13.-  Sheai>-stress  distribution 


Figure  Hi. — Distribution  of 


Figure  1.5.  - Distribution  of  u1  near  wall, 
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Figure  21.-  Distribution  Oi  Microscales. 
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Figure  25*-  Turbulent -energy  balance  near  wall 
Dashed  iine3  are  estimated. 


Figure  26.*  Chance?.  s:id  boundary-layer  mean  velocities  compared  on  basis  of 

velocity^defect  law. 
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Figure  }2 A t.vo- dimensionality  test  of  dynamic  pressure. 


. uc 


w 

> 

u * 

degrees  1 
0 

- 1 
- ? 


I 

i 


X - meter 
Claw  probe 


1.0 


Distribution  of  W relative  to  U at  station  midplane. 


- 

o 

c 

o 

o 

U~ 

V* 

0 1.0 


•Distribution  of  f^-JL.  at  station  5>  midplane. 
U 


Figure  33.-  Quantities  to  establish  test  con-HtLons 
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Figure  37.-  Hot-wire  heat  loss  curves 
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Figure  38.-  Experimental  verification  of  King's  expression  for  A •=  AtT^) 
for  a normal  and  inclined  platinum  wire. 
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Figure  3 Hot-wire  ^ correction  ba3ed  upon  equation  (22) 
Parameter  is  V “ 30°  F. 
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Figure  UO.-  Diagrams  of  hot-wire  orientation  iw  otre 
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Constant  $ calibration 
Figure  hi,-  X-meter  calibration  curves. 
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